Observable implications of geometrical and dynamical aspects of freeze-out in heavy 

ion collisions 
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Using an analytical parameterization of hadronic freeze-out in relativistic lieavy ion collisions, 
we present a detailed study of the connections between features of the freeze-out configuration 
and physical observables. We focus especially on anisotropic freeze-out configurations (expected 
in general for collisions at finite impact parameter), azimuthally-sensitive HBT interferometry, and 
final-state interactions between non-identical particles. Model calculations are compared with data 
taken in the first year of running at RHIC; while not perfect, good agreement is found, raising 
the hope that a consistent understanding of the full freeze-out scenario at RHIC is possible, an 
important first step towards understanding the physics of the system prior to freeze-out. 
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I. INTRODUCTION 

The first data from collisions between heavy nuclei at 
the Relativistic Heavy Ion Collider (RHIC) have gener- 
ated intense theoretical efforts to understand the hot, 
dense matter generated in the early stage of the colli- 
sion . Testing these theoretical ideas relies on compar- 
ison to experimental observables. Leptonic or electro- 
magnetic observables are believed to probe directly 
the early, dense stage of the collision. Most of the early 
data from RHIC, however, have been on hadronic ob- 
servables. Measurements of hadrons at high transverse 
momentum {px) !3| have generated much excitement, as 
they may provide useful probes of the dense medium pro- 
duced at RHIC 0. However, the medium itself decays 
largely into the soft (low-p^) hadronic sector. 

Soft hadronic observables measure directly the final 
"freeze-out" stage of the collision, when hadrons de- 
couple from the bulk and free-stream to the detectors. 
Freeze-out may correspond to a complex configuration in 
the combined coordinate-momentum space, with collec- 
tive components (often called "flow") generating space- 
momentum correlations, as well as geometrical and dy- 
namical (flow) anisotropics. A detailed experimental- 
driven understanding of the freeze-out configuration is 
the crucial first step in understanding the system's prior 
evolution and the physics of hot colored matter. 

In this paper, we explore in detail an analytic pa- 
rameterization of the freeze-out configuration, which in- 
cludes non-trivial correlations between coordinate- and 
momentum-space variables. We discuss the connections 
between the physical parameters of the model and ob- 
servable quantities. If the model, with correct choice 
of physical model parameters, can adequately repro- 
duce several independent measured quantities, then it 
might be claimed that this "crucial first step," mentioned 
above, has been performed. 

A consistent reproduction of all low-px observations 
at RHIC is not achieved in most physical models which 



aim to describe the evolution of the collision. In partic- 
ular, it is difficult to reproduce momentum-space mea- 
surements while simultaneously describing the freeze- 
out coordinate-space distribution probed by two-particle 
intensity interferometry measurements (also known 
as Hanburry-Brown-Twiss or HBT measurements). 
Hadronic cascade models predict a too weak momentum 
azimuthal anisotropy and too large source sizes 0. Hy- 
drodynamic transport models describe successfully trans- 
verse mass spectra and elliptic flow but fail at describing 
pion source radii 0; some hydrodynamic models have 
successfully reproduced pion source radii but only 
with different model parameters than those used to re- 
produce spectra and elliptic flow pd| . Similarly, sophisti- 
cated hybrid transport models (e.g. AMPT jl^) require 
different model paramete rs [ 131 to reproduce data on el- 
liptic flow |l3 and HBT [16|. Good reproduction of ob- 
served values has been acheived in models which adjust 
parameters to fit data within a given freeze-out scenario, 
such as in the Buda-Lund hydro approach T7| . The work 
presented here falls into this latter category. 

The parameterization used in this paper ("blast- wave 
parameterization") is similar in form to the freeze- 
out configuration obtained from hydrodynamic calcula- 
tions but we treat the physical parameters of the 
configuration (e.g. temperature) as free parameters. Our 
main goal is simply to quantify the driving physical pa- 
rameters of freeze-out at RHIC, and the dependence of 
observables on these parameters. 

Further motivation for exploring freeze-out configura- 
tions of the type discussed here, is that they implicitly as- 
sume a "bulk" system which may be described by global 
parameters (temperature, flow strength, etc). Discus- 
sions of a "new phase of matter" and its "Equation of 
State" are only sensible if indeed such assumptions hold. 
Comparison of blast-wave calculations with several in- 
dependent measurements, then, is a crucial consistency 
check of these assumptions (though, of course, a success- 
ful comparison still would not constitute a proof of their 
validity) . 
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In transport models, whether the constituents are 
hadrons Fg. l20l | , partons , or fluid elements 0, 

I2E 1^ l25||. if they re-interact substantially, pressure gra- 
dients are generated, leading to collective velocity fields 
( "flow" ) , pushing the matter away from the hot center of 
the collision and into the surrounding vacuum. Evidence 
of collective flow, generated by final state re-interaction 
of collision products, has been based largely on inter- 
pretations of transverse mass spectra and transverse mo- 
mentum azimuthal anisotropy Q . However, this scenario 
has been challenged by new measurements of p-Au col- 
lisions (2^ and new theoretical interpretations In- 
deed, so-called initial state effects such as random walk 
of the incoming nucleons |2^ or Color Glass Condensate 
phenomena |27| may offer an alternative explanation of 
the measured spectra and anisotropics in transverse mo- 
mentum. This ambiguity apparently threatens the con- 
cept that a bulk system has been created at all. However, 
it is important to recall that collective expansion, if it ex- 
ists, would manifest itself not only in momentum-space 
observables, but would also generate space-momentum 
correlations, which can be measured via two-particle cor- 
relations. 

The possible validity of any scenario may only be 
claimed if a single set of model parameters allows a suc- 
cessful description of all measured observables. Here, 
we study, in the context of a bulk collective flow sce- 
nario, transverse momentum spectra, momentum-space 
anisotropy ("elliptic flow"), HBT interferometry, and 
correlations between non-identical particles. 

Similar studies have been reported previously jlE HE 
li^lsSE^I- aspects in our study include: consider- 

ation of a more general (azimuthally-anisotropic) freeze- 
out configuration, applicable to non-zero impact param- 
eters; model studies of azimuthally-sensitive HBT inter- 
ferometry and correlations between non-identical parti- 
cles; and a multi-observable global fit to several pieces of 
published RHIC data. 

This paper is organized as follows: In Section II, we 
describe the blast-wave parameterization. In Section HI, 
we investigate in detail the sensitivity of several observ- 
ables {pt spectra, elliptic flow, pion HBT radii, and 
average space-time separation between different particle 
types) on the physical parameters of the blast-wave pa- 
rameterization. In Section IV, we perform fits to pub- 
lished data measured at RHIC for Au-|-Au collisions at 
y/sNN = 130 GeV, and, based on these fits, describe how 
as-yet unpublished analyses (azimuthally-sensitive HBT 
interferometry, and correlations between non-identical 
particles) are expected to look. The reader primarily 
interested in the quality of the fit to the data and the 
resulting parameters may want to skip past the details 
of Section HI. In Section V, we summarize and conclude 
on the relevance of the blast-wave parameterization at 
RHIC. 



II. THE BLAST-WAVE PARAMETERIZATION 

A. Geneology and Motivation 

More than a quarter of a century ago, Westfall et al [s^ 
introduced the nuclear fireball model to explain midra- 
pidity proton spectra. The idea was that the overlapping 
nucleons of the target and projectile combined to create 
a hot source with velocity between that of the target and 
projectile. Protons emitted from this source were ex- 
pected to be emitted isotropically with a thermal energy 
distribution. 

Soon thereafter, Bondorf, Garpman, and Zimanyi [3^ 
derived a non-relativistic expression for the energy spec- 
tra of particles emitted from a thermal exploding source. 
The radial fiow in their (spherical) source results in en- 
ergy spectra increasingly different than those from a 
purely thermal (non-flowing) source jas the particle mass 
increases. Siemens and Rasmussen [23 then generalized 
the formula with relativistic kinematics, further simpli- 
fying by assuming a single expanding radial shell. 

While a spherically-expanding source may be expected 
to approximate the flreball created in lower-energy colli- 
sions, at higher energies stronger longitudinal flow may 
lead to a cylindrical geometry. A decade ago, Schned- 
ermann et al [s^ introduced a simple functional form 
for the phasespace density at kinetic freezeout, which 
approximated hydrodynamical results assuming boost- 
invariant longitudinal flow [s^l, and successfully used it 
to fit Pt spectra with only two parameters: a kinetic 
temperature, and a radial fiow strength. The coordinate 
space geometry was an infinitely long solid cylinder (and 
so should approximate the situation for 6 = collisions 
at midrapidity); the transverse radial flow strength nec- 
essarily vanished along the central axis, and is assumed 
maximum at the radial edge. Most hydrodynamic calcu- 
lations yield a transverse rapidity flow field linear in the 
radial coordinate |^| . 

Huovinen et al Il8l | generalized this parameterization 
to account for the transversely anisotropic fiow field 
which arises in non-central collisions, and which gener- 
ates an elli ptic fiow signal similar to that seen in mea- 
surements |3q . This added one more parameter- the 
difference between the flow strength in, and out of, the 
reaction plane. The spatial geometry remained cylindri- 
cal, though it was assumed to be a cylindrical shell, not 
a solid cylinder. 

The measured elliptical flow systematics as a function 
of Pt and mass are fairly well-fit with the Huovinen pa- 
rameterization jssj . However, better fits were achieved 
when the STAR Collaboration generalized the model 
even further, adding a fourth parameter designed to ac- 
count for the anisotropic shape of the source in coordi- 
nate space [3^. A shell geometry was still assumed. 

To calculate the spatial homogeneity kngths probed by 
two-particle correlation measurements ^ , we must revert 
from the unrealistic shell geometry to a solid emission 
region (infinite series of elliptical shells). Furthermore, 
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additional parameters corresponding to the source size, 
emission time, and emission duration must be included, 
increasing the number of parameters 01 ■ A similar gen- 
eralization has been studied by Wiedemann "i?] . Finally, 
in this paper, we explore the effects of a "hard-edge" 
versus a smooth spatial density profile; similar studies 
have recently been done by Tomasik et al j45| and Peitz- 
mann [3Ql] for the more restricted case of a transversely- 
isotropic source. This brings to eight the total number 
of parameters which we study. 

Although the blast-wave functional form was moti- 
vated by its similarity to the freeze-out configuration of 
a real dynamical model (i.e. hydrodynamical solutions), 
it is not necessarily true that the hydrodynamical freeze- 
out configuration corresponds to the parameter set that 
best describes the data. In this sense, the blast-wave 
model presented here remains only a parameterization. 
With eight freely-tunable parameters, it is clearly a toy 
model with little predictive power. However, the goal is 
to see whether a consistent description of the data from 
the soft sector at RHIC is possible within a simple boost- 
invariant model with transverse collective flow. If this 
turns out to be the case, then it is worthwhile consid- 
ering that the parameter values indeed characterize the 
size, shape, timescales, temperature, and flow strengths 
of the freeze-out configuration. A consistent parameter- 
ization in terms of such physical quantities represents a 
true step forward, and provides valuable feedback to the- 
orists constructing physical models of the collision. 



B. Parameters and Quantities in the blast-wave 



The eight parameters of the blast-wave parameteriza- 
tion described in this paper are T, po: P2, Ry, Rx, as, tq, 
and At; their physical meaning is given below. 

The freeze-out distribution is infinite in the beam (z) 
direction, and elliptical in the transverse (x — y) plane. 
(The X — z plane is the reaction plane.) The transverse 
shape is controlled by the radii Ry and R^, and the spa- 
tial weighting of source elements is given by 



r!(r,0,) = f7(f) 



1 



1 + e(''-i)/'^= 

where a fixed value of the "normalized elliptical radius 
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FIG. 1: (Color online) The source weighting function f2, as a 
function of the normalized elliptical radius r for several values 
of the surface diffuseness parameter as. 



bution. In particular, as we discuss especially in Sec- 
tions 1111 CI and nil Ul non-zero collective fiow induces 
space-momentum correlations which dominate the spa- 
tial source density distributions. Only for a system with- 
out flow {po = /52 = 0, see below), the source distribution 
is given by fl, so that, e.g., for = 0, there is a uniform 
density of sources (^^ = const) inside the ellipse de- 
fined by Ry and R^, and no sources outside. 

The momentum spectrum of particles emitted from a 
source element at {x, y, z) is given by a fixed temperature 
T describing the thermal kinetic motion, boosted by a 
transverse rapidity p{x,y). This is common in models of 
this type. However, unlike transversely isotropic parame- 
terizations, the azimuthal direction of the boost (denoted 
(/){,) is not necessarily identical to the spatial azimuthal 
angle 05 — tan~^(y/x). Instead, in our model, the boost 
is perpendicular to the elliptical sub-shell on which the 
source element is found; see Figure 13 We believe this 
to be a more natural extension of an "outward" boost 
for non-isotropic source distributions than that used by 
Heinz and Wong (sj , who used an anisotropic shape but 
always assumed radial boost direction {(f>b = (f>s)- It may 
be shown that for our model 



corresponds to a given elliptical sub-shell within the solid 
volume of the freeze-out distribution. 

The parameter as corresponds to a surface diffuseness 
of the emission source. As shown in Figure ^ a- hard 
edge ("box profile") may be assumed by setting = 0, 
while the density profile approximates a Gaussian shape 
for Og w 0.3. 

It should be noted that the weighting function 
VL{r,(j)s) is not, in general, the source density distri- 



tan(0,) = ( ^ 



tan(0f,). 



(3) 



Hydrodynamical calculations for central collisions (i.e. 
azimuthally isotropic freezeout distribution) suggest that 
the flow rapidity boost depends linearly on the freeze-out 
radius 24!|. We assume a similar scenario, but in our 
more generalized parameterization, the boost strength 
depends linearly on the normalized elliptical radius f de- 
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isotropic boost (p2 = 0) [46 
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FIG. 2: Schematic illustration of an elliptical sub-shell of the 
source. Here, the source is extended out of the reaction plane 
[Ry > Rx)- Arrows represent the direction and magnitude of 
the flow boost. In this example, p2 > (see Equation llj . 



fined in Equation|21 Thus, in the absence of an azimuthal 
dependence of the flow (to be introduced shortly), all 
source elements on the outer edge of the source boost 
with the same (maximum) transverse rapidity po in an 
"outward" direction. 

In non-central collisions, the strength of the flow boost 
itself may depend on azimuthal angle, as suggested by 
Huovinen et al . As those authors did, we incorporate 
this via a parameter p2 , which characterizes the strength 
of the second-order oscillation of the transverse rapidity 
as a function of (jjh. Hence, the flow rapidity is given by 



p(r, (jjs) =f{po + p2 cos(20h)) , 



(4) 



We note that source anisotropy enters into our parame- 
terization in two independent ways, and each contributes 
to, e.g., elliptic flow. Setting p2 > means the boost is 
stronger in-plane than out-of-plane, contributing to pos- 
itive elliptic flow. However, even if p2 = (but po 7^ 0), 
setting Ry > still generates positive elliptic flow, since 
this means there are more sources emitting in-plane than 
out-of-plane (see Figure The STAR Collaboration 
found that both types of anisotropy were required to fit 
their elliptic flow data [s^. In generalizing the circu- 
lar transverse geometry parameterization of Huovinen et 
al (in which (j)b = 0s), they added a parameter S2 
and weighted source elements with a given (j)), as 

dN dN ,^ ^ , 

= :7:r ^ (i + 2s2cos(20s)) (5) 

Thus, a positive value of S2 = (cos(206)) corresponded to 
more source elements emitting in-plane, similar to setting 
in our parameterization. 
To facilitate comparison of fits with the STAR model 
and with ours, we relate the S2 of STAR to the geomet- 
ric anisotropy of our parameterization. In the case of 



Ry > R 



S2 



(cos(20b)) = - 



(6) 



If P2 7^ 0, anisotropics in the space-momentum correla- 
tions lead to a significantly more complicated expression. 

Finally, since our model is based on a longitudinally 
boost-invariant assumption, it is sensible that the freeze- 
out occurs with a given distribution in longitudinal 
proper time r = ^Jt"^ — . We assume a Gaussian distri- 
bution peaked at tq and with a width At 



dN 
d7 



exp 



2At2 



(7) 



We note that although the source emits particles over 
a finite duration in proper time r, we assume that none 
of the source parameters changes with r. This is obvi- 
ously an oversimplification valid only for small At; with 
time, one may expect the flow field to evolve (increase or 
decrease), and it is natural to expect the transverse sizes 
Rx and Ry to change (grow or fall) with time. However, 
calculation of the time dependence of these parameters 
requires a true dynamical model and is outside the scope 
and spirit of the present work. 



C. The Emission Function 

Our emission function is essentially a generalization of 
azimuthally- iso trop ic emission functions used by previous 
authors [s^, EH EE E3 ' ^"^^ ^iS^B we follow closely js^ : 

S{x,K) = niT cosh{ri ~ Y)Q{r,4's)e 2Ax^ x 
1 



^K-u/T ± I 

= Tot cosh(77 — F)ri(r, 0s)e sat^' x (8) 
00 

y^(zpl)"+ie~"-'^'"/'^. 

n=l 

Where the upper (lower) sign is for fermions (bosons). 
Often, only the first term in the sum in Equation [S] 
is used, resulting in a Boltzmann distribution for all 
particles. Below, we show that there is a small 
change to observables when truncating after the sec- 
ond term, and negligible effect when including further 
terms. The Boltzmann factor exp(— if • u/T) arises from 
our assumption of local thermal equilibrium within a 
source moving with four- velocity u^(a;). We assume 
longitudinal boost invariance by setting the longitudi- 
nal flow velocity vj^ — z/t (z = Tsinhry and t — 
Tcosh?/), so that the longitudinal flow rapidity ryaow = 
^ In [(1 + vl)/{1 — vl)] is identical [s^l to the space-time 
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rapidity r] = ■^ln[{t + z)/ [t — z)]. Thus, in cylindrical 
coordinates 

Ufi{x) — {coshrj cosh p{r, (l)s), sinh p(r,(j)s) cos (j)b, 

sinh p{r, (j)s) sin (j)b, sinh?7C0shp(r, 4>s)) (9) 



and 



= {niT cosh Y, px cos 4>p,pT sin 4>p, rtiT sinh Y) , (10) 

where the transverse momentuni (pt), transverse mass 
(tot), rapidity (y), and azimuthal angle {(j^p) refer to 
the momentum of the emitted particle, not the source 



element. (Note that three azimuthal angles- 
(j)p- are relevant to this discussion.) Thus 



and 



K^u^ = niT cosh p{r, (ps) cosh{ri — F) — 

Pt sinh p{r, (jis) cos{(j)b - (j)p), (11) 

and the emission function (EquationlSJ may be rewritten 
as 



S{x,K) = S{r,(f)s,T,r]) 



rriT cosh{r] — Y)il{r,(j)s)e 2Ar^ x (12) 



^(Tl)"+^e 



na cos(05 — 0p)^— n/3 cosh (77— y) 



n=l 

where we define 



Pt 

" = ^sinhp(r, 0s) 
(3= — coshp(r, 



(13) 
(14) 



Exploiting the boost invariance and infinite longitudi- 
nal extension of our source, and focusing on observables 
at midrapidity and using the longitudinally co-moving 
system (LCMS) for HBT measurements, we may simplify 
Equation El by setting Y — 0. 



where the r and rj integrals are denoted 



Ho 
Hi 

H2 



, ,-11 -(^-^0) 



= \/27r • Ar • To 

= ■ At ■ (At^ + T^) 

= V2^ ■ At-tq (SAt^ + T^] 



(16) 



and 



GoAx,K) 

Go,2(x,X) 
G2fi{^,K) 



/OC 
-00 



= 2 



2Ki(/3) 
Ki(/3) 



Ko(/3) 



K2(/3) 



(17) 



/3 

K2(/?) 

/3 ' 



+ Ki(/3) 



/3 was defined in H14|l . and Kn are the modified Bessel 
functions. For the above, we define the notation 

°o r /•2tt i-oc 

{B'}^,,(i^)^^ (Tl)"+W dct>s rdr[ 

n=l ^ -^0 

GjAx,nK)B'{x,K) X e''°"'°''^''"'"l'''^n{r,(j>s) J (18) 



for the remaining integrals, which we perform numeri- 
cally. (Note that Gij{x,nK) retains dependence on r 
and (j)s due to its dependence on f3, as defined in Equa- 
tion ll4l and so cannot move outside the integrals in Equa- 
tion HI ) 



D. Calculating Observables 

All observables which we will calculate are related to 
integrals of the emission function H12|l over phasespace 
d'^x = dxdydzdt = rdrdrjrdrdcfis , weighted with some 
quantity B(x, K). In all cases, the integrals over r and 77 
may be done analytically, though the result depends on 
whether B{x,K) itself depends on t and ry. 

In particular, if B{x,K) = 

B'{r,(j)s,K)T^ smh-' 7] cosh'^ T] then the integrals of 
interest are 

^<^^ /o°° ^'^^ iZo '^^ iZo ^'^'^^ K)B{x, K) 

= mrH, ■ {B'}^ , (K) (15) 



III. CALCULATION OF HADRONIC 
OBSERVABLES 

In this Section, we discuss how hadronic observables 
are calculated from the parameterized source and illus- 
trate the sensitivity of these observables to the various 
parameters presented in Section III Bl 

With several observables depending on several param- 
eters, it is not feasible to explore the entire numerical 
parameter space. Instead, we anticipate the results of 
the next Section, in which we fit our model to existing 
data, and vary the parameters by "reasonable" amounts 
about values similar to those which fit the data. The de- 
fault parameter values used in several of the calculations 
in this Section are listed in Tabled 
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parameter round source non-round source 










0.05 




12.04 




11 


Hy (tm) 


iz.U4 




1 Q 


T (GeV) 




0.1 




Po 




0.9 




TO (fm/c) 




9 




At (fm/c) 




2 




as 










TABLE I: Default parameter values for most calculations in 
Section irm Note that p2, Rx and Ry default values depend on 
whether we are discussing an azimuthally isotropic ( "round" ) 
or anisotropic ("non-round") source. One might expect such 
sources from central and peripheral collisions, respectively. 



A. Pt spectra 




1.5 : 

Pt (GeV/c) 



In the notation of Section III Dl the (azimuthally- 
integrated) pt spectrum is calculated as 

= [ d(t)p [ d'^xS(x,K) 
PTdpT J J 

^ niT f d0p{l}o o(^)- (19) 



FIG. 3: Transverse momentum spectra for protons (upper 
curves) and pions (lower curves), as calculated by Equa- 
tion 1191 for several values of A'^, the maximum value of n 
taken in the summation of Equation |H| see text for details. 
Parameters values correspond to the "non-round" source of 
Table |I] All spectra are arbitrarily normalized to unity at 
Pt = 0. 



In this paper, we focus only on the shapes, not the nor- 
malizations, of the spectra. 

We note that spectra calculated in the blast-wave 
model scale neither with uit nor pT, as both quantities 
enter the expression through a and /3 (Equations 1)1 3|l 
and 1)14(1 ). This breaking of mr-scaling is a well-known 
consequence of finite transverse fiow |3y| (p ^ in our 
model) . 

According to Equation ^1 rriT spectra calculated in 
the blast-wave model are insensitive to the time param- 
eters To and At. The spectral shapes are furthermore 
insensitive to the spatial scale (i.e. Ry) of the source, 
though, as we see below, there is some small sensitivity 
to the spatial shape (i.e. ^). 

First, we study the importance of using quantum (as 
opposed to classical) statistics in the source function. 
Figure 121 shows pr spectra for pions and protons, treated 
as bosons and fermions, respectively. Model parameters 
were set to the "non-round" values listed in Table The 
sum in Equation |S1 (and Equations E| and I18|l is over 
n = 1 . . . N; curves are shown for N — 1,2, 3, 4. For pa- 
rameter values in the range we study here, proton spec- 
tra are essentially independent of N. For iV > 1, the 
pion spectra are likewise robust against the value of N, 
though in the classical limit {N — 1), there is relatively 
lower yield at low pr- (Note that all spectra are arbitrar- 
ily normalized to unity at = 0.) Calculations below 
use the truncation N = 2. 



1. Spectra from Central Collisions 

Focusing first on central collisions (so that the flow 
anisotropy parameter p2 = and Rx = Ry), then, we 
need only consider the spectra sensitivity to the temper- 
ature and radial flow parameters po and T, and to the 
surface diffusion a^. 

Fixing the transverse spatial density distribution to a 
box profile (og = 0), the evolution of spectral shapes 
for pions and protons are shown in Figures 01 and El as 
the temperature and radial flow parameter, respectively, 
are varied about nominal values of T = 100 MeV and 
Po = 0.9. As has been noted previously [s^l, at low 
Pt, temperature variations affect the lighter pions more 
strongly, while variations in the collective flow boost pro- 
duce a stronger effect on the heavier particles. 

Next, we consider the effect of a finite surface diffuse- 
ness parameter (os ^ 0), i.e., using the smoother spa- 
tial density distributions of Figure ^ Since we assume 
a transverse flow proflle which increases linearly with ra- 
dius (c.f. Equation 01), one trivial effect is that, for a 
fixed pq, increasing will produce a larger average fiow 
boost (p). The effect of increasing transverse flow was 
already explored directly in Figure |S1 so we avoid this 
trivial effect here, and explore the effect of varying a^, 
while keeping {p) constant [30l |. 

The average transverse flow boost (p) = F{as) ■ po 
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FIG. 4: Transverse momentum spectra for protons (upper 
curves) and pious (lower curves), as calculated by Equa- 
tion 1191 for several values of the temperature parameter T. 
Other parameters follow the "round" source defaults of Ta- 
ble AH spectra are arbitrarily normalized to unity at 
Pt ~ 0. 




Pt (GeV/c) 

FIG. 5: Transverse momentum spectra for protons (upper 
curves) and pious (lower curves), as calculated by Equa- 
tion 1191 for several values of the radial flow parameter po • 
Other parameters follow the "round" source defaults of Ta- 
ble AH spectra are arbitrarily normalized to unity at 
Pt = 0. 



FIG. 6: The geometric constant of proportionality F between 
the average transverse flow boost (p) and the blast-wave pa- 
rameter po, as a function of the surface diffuseness as. 




Pt (GeV/c) 

FIG. 7: Transverse momentum spectra for protons (up- 
per curves) and pions (lower curves), as calculated by Equa- 
tion^! for several values of the surface diffuseness parameter 
Os. The radial flow strength po is co- varied; see text for de- 
tails. Other parameters follow the "round" source defaults 
of Table Q] All spectra are arbitrarily normalized to unity at 
Pt = 0. 

where the geometric proportionaUty constant 



F{a,) = l+exp((.-l)/a O (20) 

Jo l+oxp((x-l)/a.,) 

is independent of po or R — — Ry. For the box profile, 
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F{as = 0) = |. Figure shows this geometric factor as 
a function of the surface diffuseness. 

Figure [3 shows the pion and proton spectra for var- 
ious values of Us- The radial flow strength po was co- 
varied with fls so that the average transverse flow boost 
was (p) = 0.6. To a first approximation, the spectral 
shapes depend only on the temperature T and the aver- 
age transverse flow boost (p). The residual dependence 
on the surface diffuseness parameter Og arises from the 
fact that while the average boost rapidity has been held 
constant, the spread of boost rapidities increases with in- 
creasing Qs |30j. Thus, we observe qualitatively similar 
variations in the spectral shapes when as increases (Fig- 
ure Ej), as when T increases (Figure 0J. The variations 
are not quantitatively identical since in the present case, 
the velocity spread is not thermal, and the particle ve- 
locity spread evolves differently with mass, depending on 
whether it arises from a boost spread or a thermal spread. 



2. Dependence of spectral shapes on source anisotropy 

Azimuthally-integrated pr spectra are often presented 
as a function of event centrality. For 6^0 collisions, the 
emitting source may have anisotropic structure {Rx ^ Ry 
and p2 7^ in the present model). Thus, it is interesting 
to explore possible effects of these anisotropics. 

For an azimuthally isotropic flow field (p2 = 0), the 
spectral shapes are insensitive to spatial anisotropics in 
the source (i.e. Rx Ry)- This is because the spectral 
shapes are determined by the distribution of boost veloc- 
ities |3Cll |. which is unchanged by a shape change in our 
parameterization, if p2 — 0. 

For an azimuthally-symmetric spatial source {Rx = 
Ry), a very small variation in the (/)p-integrated spectral 
shapes is observed when p2 is changed from a value of 
0.0 to 0.15, as seen in Figure |S1 This, again, is due to 
the slightly increased spread in boost velocities; source 
elements emitting in-plane boost a bit more, and out-of- 
plane a bit less. This effect becomes stronger in the pres- 
ence of an out-of-plane spatial anisotropy {Ry/ Rx > 1) 
as shown in Figure El In any case, the effects of "rea- 
sonable" source anisotropy on the shapes of azimuthally- 
integrated spectra are very small. 

Thus, we conclude that azimuthally-integrated pt 
spectra are largely insensitive to "reasonable" source 
anisotropics (see Section HVI for "reasonable" ranges) and 
probe mainly the thermal motion (T) and average trans- 
verse fiow boost ((p)) of the source. 



B. Elliptic flow versus mass and pt 

In the notation of Section III DI thc elliptic fiow param- 
eter V2 is calculated as 

, , Iq" d<Pp{cos{2cl)p)}^ ^{K) 
V2{pT,m) = — (21) 




1.5 2 

Pt (GeV/c) 

FIG. 8: Transverse momentum spectra for protons (upper 
curves) and pions (lower curves) , as calculated by Eauation ll9l 
for an azimuthally symmetric flow field (p2 ~ 0) and an asym- 
metric field (p2 ~ 0.15). Other parameters follow the "round" 
source defaults of Tabled] AH spectra are arbitrarily normal- 
ized to unity at pT = 0. 
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FIG. 9: Transverse momentum spectra for protons (up- 
per curves) and pions (lower curves), as calculated by Equa- 
tion^] for an azimuthally anisotropic flow field {p2 = 0.15), 
and a spatially isotropic {Ry/Rx = 1) and anisotropic 
{Ry/ Rx ~ 1.5) spatial distribution. Other parameters follow 
the "round" source defaults of Table 13 



A finite V2 arises from azimuthal anisotropics in the 
source {Rx ^ Ry and/or p2 ^ 0). As discussed below, 
however, the parameters po and T strongly affect its value 
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FIG. 10; Elliptic flow parameter «2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for several values of 
A^, the maximum value of n taken in the summation of Equa- 
tion|H| see text for details. Model parameters follow the "non- 
round" source defaults of Table IH 



FIG. 11; Elliptic flow parameter «2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the temperature parameter T. Other parameters follow the 
"non-round" source defauhs of Table |I| 



and evolution with px and mass. In the present param- 
eterization, V2 is not sensitive to the overall spatial scale 
of the source (Ry) or the time parameters tq and At. 

The V2 parameter depends non-trivially on both px 
and particle mass 0, US this section, we 

explore the evolution of pion and proton V2, as we vary 
the model parameters from nominal "non-round" values 
for non-central collisions (cf^J. 

As with the px spectra of Section IIII Al we first check 
the importance of quantum statistics. Figure IIUI shows 
V2 for pions and protons, for different values of N, where 
the sum in Equation |S1 (and Equations El and I18|l runs 
over 71 = 1 . . . A^. Again, we find only a small difference 
for the pions between A^ = 1 (classical limit) and N = 2, 
beyond which V2 is robust against further increases in A^. 
Calculations here use N = 2. 

Figure E] shows the evolution of V2 as the tempera- 
ture parameter (T) is varied. For both particle types 
shown, the increased thermal smearing in momentum 
space, as T is increased, leads to a reduced momentum- 
space anisotropy. The effect of the thermal smearing is 
greater for the lighter pions. 

Less intuitive is the evolution of V2 as the flow field 
{po or P2) is varied. In Figures [T^ and ICT the average 
transverse flow parameter po is varied for an azimuthally 
isotropic {Ry/Rx = 1) and anisotropic {Ry/Rx = 13/11) 
shape, respectively. For the isotropic spatial distribution, 
we find that V2 decreases as po increases, for all px and 
for both particle types. This is due to the decreasing 
relative amplitude of the oscillation in the flow field. Not 



■0.12 
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Pt (GeV/c) 



FIG. 12: Elliptic flow parameter V2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the transverse flow parameter po. The source spatial distribu- 
tion is assumed azimuthally-anisotropic {Ry — R^), but other 
parameters follow the "non-round" source defaults of Table HI 



surprisingly, the effect is larger for the heavier protons. 

Indeed, it has been pointed out 'i'S, 'ilj that high 
radial flow (large po) can lead to negative values of V2 
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FIG. 13; Elliptic flow parameter «2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the transverse flow parameter po- Other parameters follow 
the "non-round" source defaults of Table El 



FIG. 15: Elliptic flow parameter «2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the modulation in the transverse flow p2. Other parameters 
follow the "non-round" source defauhs of Table |T| 
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FIG. 14: Elliptic flow parameter «2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the modulation in the transverse flow p2. Other parameters 
follow the "round" (note: Ry = R^^) source defaults of Table|I] 



for heavy particles; this is clearly true for the protons 
in Figure El This negative V2 reflects the depletion of 
the \ow-pT particle yield when most source elements are 
highly boosted transversely; this is also the origin of the 



"shoulder arm" in the pr spectra (c.f. Section fill Ap . 
This depletion is larger in-plane (since the boost is higher, 
i.e. p2 > 0), leading to negative V2- At higher pT, and/or 
lower flow strengths, the competing effect is dominant: 
the larger in-plane boost leads to more particles emitted 
in-plane, and V2 > 0. 

For a slightly anisotropic shape {Ry/Rx = 13/11; Fig- 
ure the V2 parameter increases significantly for both 
protons and pions, due to the larger number of source 
elements boosted in-plane (c.f. FigureEJ. The finite spa- 
tial asymmetry leads to an effect that tends to oppose 
the reduction of V2 with increasing poj discussed above, 
and, for pions at low px, even reverses it in this case. 
Hence, at low pT, increasing pQ increases (decreases) V2 
for pions (protons). 

We observe similar trends when po is held fixed, and 
P2 is varied. Figure FPil shows V2 evolution for a spatially 
aziniuthally-symmetric source {Ry/Rx — 1) and Fig- 
ureElfor a slightly asymmetric source {Ry/Rx = 13/11). 

As mentioned above, a finite value of V2 may be ob- 
tained for an azimuthally- symmetric flow field, if the spa- 
tial shape is asymmetric. Figure 1161 shows the proton 
and pion V2 for p2 = and po = 0.9, for various values 
of Ry/Rx- The larger number of sources emitting in- 
plane results in more particles measured in-plane (thus 
V2 > 0). This momentum-space asymmetry saturates 
at Pt ^ 1.0 GeV/c (for this set of parameters) and is 
relatively independent of particle mass. (We note that 
the pT-scale in Figure ^1 is larger than the other fig- 
ures, in order to show the V2 saturation effect.) This 
effect is similar to those previously discussed by Hou- 
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FIG. 16; Elliptic flow parameter «2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the spatial anisotropy (Ry/Rx). Elliptic flow parameter V2 for 
pions (upper curves) and protons (lower curves), as a func- 
tion of transverse momentum, as calculated by Equation 1211 
for various values of the modulation in the transverse flow p2. 
Other parameters follow the "round" (note: p2 = 0) source 
defauhs of Table 13 



FIG. 17; Elliptic flow parameter V2 for pions (upper curves) 
and protons (lower curves), as a function of transverse mo- 
mentum, as calculated by Equation 1211 for various values of 
the surface diffuseness parameter Os. The values of po and 
P2 are co-varied with as so that the average transverse flow 
and its average azimuthal modulation remain flxed; see text. 
Other parameters follow the "non-round" source defaults of 
Table □ 



vinen, et al, [s^ and Shuryak js^- However, we stress 
the non-zero V2 parameter does not indicate "elliptic flow 
without transverse flow," as suggested in |^ . If trans- 
verse flow is turned off (po = 0), the spectra are thermal 
and isotropic, and W2 — for all pT] space-momentum 
correlations (induced by p 7^ in our model) are re- 
quired for a finite V2. In the Houvinen model [s^ it is 
the implementation of the Cooper- Frye freeze-out proce- 
dure (which creates an infinitely opaque source) which 
produces the space-momentum correlations, effectively 
generating "flow." Shuryak implemented a parame- 
ter K which controlled the opaqueness; V2 was found to 
increase with n until the source was essentially infinitely 
opaque. 

Finally, we again explore the effect of "softening" the 
spatial source distribution r2(r, 0s). Similar to the dis- 
cussion surrounding Figure it is clear that the most 
meaningful comparison comes from keeping the average 
value of the transverse boost (and its anisotropy) con- 
stant, not keeping the parameters po and p2 constant, as 
we vary Os . Similar to the discussion of the pr spectra, 
we find in Figure [T7I that V2 is largely insensitive to the 
surface diffuseness parameter Os. 



C. HBT radii versus px and (j>p 



The previous subsections have discussed momentum- 
space observables only, even though coordinate-space 
considerations came into play indirectly. The most di- 
rect experimental probes of the space-time structure of 
the freeze-out configuration are two-particle correlation 
functions in relative momentum |^, |63 • Here, we use the 
standard "out-side-long" coordinate system of Pratt and 
Bertsch in which the long direction (i?/) is parallel to 
the beam, the side direction (i?s) is perpendicular to the 
beam and total pair momentum, and the out direction 
(i?o) is perpendicular to the long and side directions. 



For boost-invariant sources, the HBT radii i?^; and i?^; 
vanish by symmetry |4?|. (For the more general case, see 
References 

|3lillilE3i3-) Thus, we are left with 
four HBT radii, which are related to space-time variances 
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as 



Rf 



where 



-{xy) sm{2(j)p) 

+ {xy) sin(2(/)p) — 2fi^{{tx) cos(j)p+{iy) sincfip) 

{iy) cos(20p) - i ((£2)_(y2)) sm{2(j)p) 
+P±{{ix) sin (l)p-{iy) cos^p), 
(z2)_2A(fe)+A'(i'>, 



(22) 



(/(^)>(i^) 



/ d'^xf{x)S{x,K) 
J d^xS{x, K) 



(23) 



We restrict our attention to correlations calculated in 
the Longitudinally Co-Moving System (LCMS), in which 
Y = (3i = Q. In this case the last Equation 1^ simplifies 
to 

i?f = {~z^) (24) 
In the notation of Section FlI Dl and further defining 



{■^b.fc = 111 ' 

the space-time correlations of interest are 



(25) 



{it) 

{yi) 

{P) 



{y^}o.o - Mo,o 



To 



To 

(3At2 



{y}o,i - {y}o,o {l}o,i 



(26) 



We note that all quantities with space-time dimensions 
are explicitly shown in Equations and in particular, 
all dependence on the timescale parameters tq and At. 

The proper time of freeze-out is often estimated 
(e.g. ,51]) by fitting the mr-dependence of the measured 
Ri radius, to a formula motivated by Sinyukov and col- 
laborators 52, 53], and subsequently improved upon by 
others They assumed boost-invariant longitudi- 

nal fiow, but vanishing transverse flow {p = ^ P — ^-^) 
and instantaneous freeze-out in proper time (i.e. At — 
0) ; they also simplified the formalism by using Boltzmann 
statistics (i.e. using only the first term in the summation 



in Equation IHl . In this case, we find, in agreement with 
References 54 .55j . 



^ K2(mT/T) ^ 
rriT Ki(mT/T) 



(27) 



unity for ^ 



rriT 
T 



We note that the last term, which represents a correc- 
tion to the original Sinyukov formula 52, 53], approaches 
oo but remains sizable (~ 1.5) for 

2. 

In general, the emission "homogeneity region" 
(characterized by the correlation coefficients of Equa- 
tion I2t)|) and the HBT radii of Equation [23 depend on 
the pair momentum 0,^3 ■ For a boost-invariant system, 
this corresponds to dependences onp^ and (j)p. Figure lTHl 
shows projections onto the transverse {x — y) plane of the 
emission regions for pions with pT = 0.3 GeV/c and 4>p — 
0°, 135°, for an anisotropic source with strong transverse 
flow. The flow generates strong space-momentum cor- 
relations, so that particles with higher px tend to be 
emitted from the edge of the source, with (j)p ^ (ps- To- 
gether with the finite extent of the overall source, this 
implies that, spatially, the emission region is often wider 
in the direction perpendicular to the particle motion (in- 
dicated by the arrows) than along the motion; this can 
have strong implications, e.g., for the difference between 
Rl and Rl. 
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FIG. 18: (Color online) Emission probability contours, plot- 
ted on a linear scale, indicate emission zones for pions with 
PT ~ 0.3 GeV/c at (l>p ~ 0° and (j)p — 135° (indicated by 
arrows), from a blast-wave source with "non-round" default 
parameters listed in Table |l] 

While the blast-wave parameterization provides direct 
access to the homogeneity region, the radii extracted 
from two-pion interferometry are compared to second or- 
der moments calculated as shown in equation 1261 Such 
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FIG. 19: Distribution of the pion spatial distribution in the 
direction x (left), y (middle) and z (right). Plain line; pion 
momentum px = 0.25 GeV/c and py = 0. Dash line: pion 
momentum p^ = 0.5 GeV/c and py = Q. Round source pa- 
rameters of Table EJwere used, with the exception of as that 
was set to 0.01. 

comparison is strictly correct only if the homogeneity re- 
gions are Gaussian distributions. Figure ^1 shows an 
example of spatial distributions along the 3 Cartesian di- 
rections for pions with two different momenta, ~p*= (0.25 
GeV/c, 0, 0) and ~p= (0.5 GeV/c, 0, 0); in this case, the 
X—, y— and z— directions correspond to "out," "side" 
and "long," respectively The z distribution is nearly 
Gaussian; the difference between the a extracted from 
a Gaussian fit and the second order moment is on the 
order of a few percent. On the other hand, the x and y 
distributions are clearly not Gaussian. To estimate the 
level of distortion introduced by the non-Gaussian shape, 
we compare the second order moments with the Gaussian 
a extracted by fitting the peak of the spatial distributons. 
Indeed, it has been shown in [g^ that calculating cor- 
relation function from models and fitting them as exper- 
imental data yields radii that are close to the Gaussian 
a extracted by fitting the peak of the spatial distribu- 
tons. In both the x— and y— directions, we find that 
the Gaussian a are systematically larger than the sec- 
ond order moments by up to 20% (depending on the fit 
range) for pions at pr — 0.25 GeV/c. This discrepancy 
diminishes when increasing the transverse momentum or 
the particle mass. E.g. it is on the order of 10% for pi- 
ons at pt = 0.5 GeV/c or for kaons at pr = 0.25 GeV/c. 
Thus, comparing second order moments with the radii ex- 
tracted from two-pion correlation functions may involve 
significant systematic errors at low transverse mass. This 
issue may be overcome by calculating correlation func- 
tions from the Blast Wave space-time distributions or by 
relying on imaging method to extract space-time distri- 
butions from the data !6lj. Applying these methods is 
beyond the scope of this paper and we will thus carry 
on keeping in mind that systematic errors are associated 
with comparing HBT radii with second order moment at 
low transverse momentum. 

On the other hand, resonance decay may introduce a 
core-halo pattern in the distribution of particle space- 



time emission points [6^. Indeed, some resonances de- 
cay sufficiently far away from the bulk of the system that 
their decay products emerge beyond the main homogene- 
ity region. Such particles form a halo around the core of 
the source. However, to affect the extracted radii, the 
resonance lifetime needs to be short enough for the cor- 
relation to take place at a relative momentum accesible 
to the experiment. This lifetime is typically between 10 
and 100 fm. The effect of resonance feed-down on radii 
measured by two-pion interferometry has been studied 
within a early version of the blast wave parameteriza- 
tion (6^ . It was found that the uj meson leads to a halo 
effect, its cr being 23.4 fm. The other resonances have 
lifetimes that are either too short (e.g. p, A) or too long 
(e.g. ?7, K^). At RHIC energy, thermal models [B^I show 
that about 10% of the pions come from w, which means 
that this effect should be rather small. The effect of the 
very long-lived (ct > 100 fm) resonances is usually as- 
sumed to reduce the so-called A parameter |65^ | , which we 
do not discuss here. 



0.6- 




p^ (GeV/c) (rod) 

FIG. 20: (Color online) Space-time correlation {xt) as a func- 
tion of for (j>p — (left panel), and as a function of 0^ for 
Pt ~ 0.3 GeV/c. "Non-round" source parameters of Table|I] 
were used. 

Also important are temporal effects, including space- 
time correlations. The average x — t correlation, quanti- 
fied by (xt), for the same anisotropic source, is plotted in 
Figure 1201 as a function oi px and cf>p. In the left panel, 
(j>p — 0°, and thus x is the "out" direction. We first note 
that the correlation is negative, tending to increase 
for finite /3± (see Eauation l22|l : large and negative "out"-t 
correlations at freezeout in some hydrodynamical mod- 
els have been a significant component of predicted 
large R^/R^ ratios. As required by symmetry (481, (xt) 
displays an odd-order cosine dependence on (j)p; clearly 
the first-order component dominates here, as shown in 
the right panel of Figure This cos(0p)-dependence is 
driven mostly by the cos{4>s) dependence of the x coor- 
dinate, coupled with the radial flow which tends to make 
0s ~ 'Pp (c-f. Figure [TH|l . However, we also note from the 
figure that the scale of these correlations is < 1 fm^/c, 
which, as we shall see, is much smaller than the typical 
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scale of R^; hence, space-time correlations do not domi- 
nate in the blast-wave. 

More important are the correlation coefficients (iy), 
(i^) and (y^), and their p-dependence. (xy) quantifies 
the transverse "tilt" of the emission zone with respect 
to the reaction plane. As required by symmetry |48| . 
this tilt vanishes at (/>p = 0, while, in the present case 
with Ry > Rx, it is positive for (pp — 135° (see Fig- 
ure lTH)) . That {x^) and (y^) also depend on (f>p is likewise 
clear from Figure [THl As we shall see, the 0p-dependence 
of {xy), (i^), and (y^), combined with the explicit (j)p- 
dcpendcnces in Equation [221 drive the oscillations in 
HBT radii which we will study. 

At this point, it is worthwhile to mention that the emis- 
sion zone ( "homogeneity region" ) and the correlation co- 
efficients ((x^), (y^), {xy), {it), {yt)) will vary with (jip 
also for an azimuthally isotropic source {Rx = Ry, P2 — 
0). However, of course, the measured HBT radii will be 
(/)p-independent. This arises due to cancellation and com- 
bination of terms in Equations 122 For i?^, the first term 
becomes (/)p-independent, and the second and third terms 
combined are 0p-independent (((i^) — (y^)) oc cos(2(/)p) 
and {xy) cx sin(2(/)p)). For R^, the story is the same for 
the ffist three terms, while the fourth term becomes (pp- 
independent as {it) oc cos((/)p) and {yt) oc sin((/)p). Mean- 
while, for the ffist and second terms cancel each 
other, as do the two components of the third term; hence 
R^g — for azimuthally-symmetric sources. We men- 
tion these points here because several of the exact can- 
cellations and combinations which hold for an isotropic 
source, continue to hold approximately even when the ge- 
ometry or flow field is anisotropic. 

We turn now to a detailed study of the observable HBT 
radius parameters, and the effects of varying blast- wave 
parameters. For a boost-invariant system, the symmetry- 
allowed tpp oscillations of the (squared) HBT radii are |43| 

RI{vtAp) = <obT) + 2^„^2 4 6,...^'«bT)cos(n(/)p), 

RI{PtAp) = Rl.oiPT)+'^En=2A,6,...KniPT)cOs{ncbp), 

RIs{ptAp) = 2E„=2A6,...^o.,JPT)sin(n<^p), (28) 

Rf{PT,(l)p) = RloiPr) +2j2n=2A,6,..RlniPT)cOs{n(l)p), 

where the HBT radii Fourier coefficients R^ „ (/^ = 
o, s, OS, I and n = 0, 2, 4, . . . ) are (/)p-independent. 

Figure |2] shows the pr and (f>p dependence of HBT 
radii calculated for a blast-wave source with a slightly 
anisotropic flow field and shape. In addition to an over- 
all decrease in the average value of the HBT radii with 
increasing pT, we observe significant oscillations in the 
transverse radii R^, Rj, R^^, and smaller oscillations in 
Rf. 

The Fourier coefficients may be calculated as 

^2 ^^^^ ^ I {RUpt, (f>p) cos(n0p)) (/i = o, s, I) 
^ \{Rl{PT,(l>p)sm{n4>p)) ifi^os) 

(29) 




FIG. 21: Squared HBT radius parameters calculated with 
Eg uations 1261 and 1221 from a blast-wave source with the "non- 
round" default parameters of Table |I] Squared radii for vari- 
ous cuts in pt are plotted versus (pp, the emission angle which 
respect to the event plane. 



In the present model, we find that oscillation amplitudes 
above 2"^^ order are very small in all cases considered 
{\Rfi,4/ Rfi,2\ ^ 0.01). Therefore, the 0p-dependence of 
the HBT radii at a given px is essentially encapsulated in 
seven numbers: the O"^- and 2"'^-order Fourier coefficients 
of Rl, and Rf, and the 2"'^-order Fourier coefficient 
of i?^^. We henceforth explore the evolution of the pr- 
dependence of these seven numbers, as model parameters 
are varied. 

Figure 1221 shows the Fourier coefficients for n = 0, 2, 
corresponding to a blast-wave source with T = 0.1 GeV, 
an isotropic flow field (po = 0.9, p2 =0), a box profile 
(fls = 0), and time parameters tq = 9 fm/c. At = 2 fm/c. 
Here, the average transverse size of the source (i?^ -I- i?^) 
was held fixed, while the shape {Ry/R^) was varied. 

The 0*'^-order Fourier coefficients (corresponding to the 
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FIG. 22: Fourier coefficients of tfie (^p-dependence of tlie 
squared HBT radii, as caicuiated by Equation 1291 0**^- and 
2"'*-order Fourier coefficients are piotted in ttie ieft and rigtit 
paneis, respectiveiy. Tlie transverse sliape of tlie source (i.e. 
Ry/Rx) was varied, whiie R^ + was lieid fixed. "Round" 
source vaiues from Tabic Q are used for tfie otlier parameters. 



radii Ro, Rs, and Ri usually measured by experimen- 
talists) are sensitive only to the average scale, not the 
shape, of the source. The average values of the trans- 
verse radii R^ n and i?^ q fall with increasing pT due to 
radial flow |6[l55l| (c.f. Figure I^T|l . At intermediate val- 
ues of PT, -Rq.o > -^5,0 due to finite timescale effects 
(cf Figures and OOjl, but at high pr, RI^q < iJ^.o 
(i.e. Ro/Rs < 1), i n q ualitative agreement with ex- 
perimental data jsj, I63. The boost-invariant longi- 
tudinal flow produces the strong decrease of i?^Q with 

PT UMMMM 

Richer detail is seen in the oscillations of the HBT 
radii, quantified by Fourier coefflcients i?^ 2 the right- 
hand panels of Figure [221 Here, the elliptical shape 
of the source is explicitly clear. The signs of the 2"*^- 
order Fourier coefficients of the transverse radii directly 
reflect the out-of-plane-extended source geometry when 
Ry > Rx- Rqs 2 has a similar geometric interpretation, 
in terms of the (/)p-evolution of the "tilt" of the homo- 



FIG. 23: (Color online) Identical data as in Figure l^ except 
the right panels show ratios of 2"'*-order and O*'' order Fourier 
coefficients of the (/>p-dependence of the squared HBT radii. 
Arrows indicate values calculated from Equation 1321 see text 
for details. 



geneity region |23|. The relatively small oscillations in 
Rf arise not directly from geometry, but instead from 
transverse flow gradients, which slightly reduce Rf [55j . 
In the present example, the transverse flow increases lin- 
early from (at the center) to pq = 0.9 (at the edge of the 
source), independent of boost angle 0&. However, when 
Ry > Rx, the flow gradient is larger for source elements 
boosted in-plane, leading to slightly greater reduction of 
Rf when cjjp — 0; hence i?^2 < 0- 

Finally, we recall that the flow field is isotropic (p2 — 0) 
and so all (/)p-dependence arises from geometry here. 
Thus, if the values of Ry and Rx are interchanged (cor- 
responding to in-plane-extended sources) , i?^ q would re- 
main unchanged, and i?^ 2 would simply change in sign 
in Figure!^ 

As discussed above, transverse flow-induced space- 
momentum correlations tend to decrease homogeneity 
lengths as px increases. When combined with other 
effects (e.g. temporal effects), non-trivial pt depen- 
dences of the HBT radii result. The pT-dependences of 
the (/)p-averaged values i?^ q have been discussed exten- 
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sively [a|55|. Meanwhile, the p^-dependences of the os- 
cillation amplitudes (i?^2) shown in the right panels of 
Figure [221 have not been explored previously and may be 
non-trivial in principle. 

It was suggested 68J that the pT-dependence of |i?^,2l 
might be driven largely by the same effects which gen- 
erate the pT-dependence of R^q, and hence the most 
efficient and direct way to study the source is to plot 
R^Q, which encode scale information, and then the ra- 
tio of 2"^^- to 0*'^-order Fourier coefficients, which would 
encode geometric and dynamic anisotropy. This is an 
excellent suggestion, though consideration must be given 
to the appropriate scaling. First, we consider the trans- 
verse radii R^, i?^ and R^^. The radii R^ and R^^ encode 
both transverse geometry and temporal information. As 
discussed above, space-time (e.g. x — t) correlations are 
small in magnitude, and furthermore affect the HBT radii 
in combinations which tend to cancel any 0p-dependence. 
Therefore, we expect i?oo, i?s_0' ^s2j -^o 2 a-^id Rls,2 to 
contain geometric contributions, while temporal contri- 
butions are significant only for R^q. In this case, the 
appropriate ratios to study are i?o^2/^s,0' ^^.2/^-1,0 ^^'^ 
R-os 2/ R^ 0- Indeed, we find numerically that these are 
the ratios least affected by the overall scale of the ho- 
mogeneity region, which varies both with px and with 
model parameter. The oscillation strength Rf2 of the 
longitudinal radius, on the other hand, is entirely due to 
implicit (pp dependences driven by space-momentum cor- 
relations; these same correlations affect -R^q. Hence, the 

appropriate ratio to study in this case is Rf2/Rfo- 

In Figure!^ we show these ratios for the same sources 
as were plotted in Figure [221 The pr-dependence of 
the ratios is significantly less than that of the oscillation 
strengths i?^ 2j ^.s anticipated, due to the fact that the 
latter is driven largely by space-momentum correlations 
reducing the spatial scale of the homogeneity region. 

Going further, we may recall that in the special case 
of vanishing space-momentum correlations (p = or 
T — 00), the transvserse radii oscillate with identical 
strengths {Rl^ 2 ~ 2 = ^-^o 2); ^.nd the in- plane and 
out-of-plane extents of the source may be directly ex- 
tracted 0147, 48, 49], as the "whole source" is viewed 
from every angle. In that special case, independent of px 

Rl,^lm + {i')) = l-{Rl + Rl) (30) 

RL2 = Rl2 = -Rl2 = \ 

= ^{Rl-Rl) (31) 

so that 

7j2 7j2 n2 ^ t>2 t>2 

Rpsq _ Rs,2 _ Ro,2 _ 1 Ry- Rx _ ^ 

R^fl R^fi R^sfi 2 Ry + R^ 2 

In the presence of flow, however, HBT radii measured 
at momentum p reflect homogeneity lengths which in 



principle may vary nontrivially both with pT and (f>p. 
While we find that non-vanishing flow violates the pt- 
independence of R^q, -R02; -^s 2 ^^^^ -R^^ 2 (^^^^^ ^hus 
Equations [20] and I31|) , Equation [221 remains remarkably 
robust. As seen in the next several Figures, the ratios 

Rh/Rlo^ Rla/R^fi Rlsa/Rlfi^ largely independent 
of pt, provide an estimate of the source ellipticity e. Ar- 
rows to the left of the panels for i?^ 2/^s oi R^ 2/Rl ^^'^ 
Ros2/Rs0 ™ Figures f23l34l indicate e/2 for the sources 
used. 

Now that we have established the quantities to be ex- 
amined in this Section, we briefly check the importance 
of using quantum, rather than classical, statistics in the 
source function of Equation[Sl Setting the parameter val- 
ues to correspond to the "non-round" source of Table ^ 
we plot in Figure the Fourier coefficients correspond- 
ing to different values of N, where the summation in 
Equation[Sl(and Eauationsll2land ll8l) is over n = 1 . . . A'^. 
Once again, we find a small difference between the curves 
for A^ = 1 and A^ > 1, while inclusion of higher terms in 
the summation have essentially no effect. Blast- wave cal- 
culations in this Section correspond to A^ = 2. 

Recently, Heinz and Kolb, in a hydrodynamic model, 
calculated HBT radii as a function of 0p for non-central 
collisions They used two different Equations of 

State and initial conditions: one ("RHIC") is appropriate 
for soft physics at RHIC energies, and has successfully 
reproduced momentum-space observables 0; the other 
("IPES") assumes extremely high initial energy densi- 
ties, perhaps appropriate for collisions at LHC energies. 

It is worthwhile to point out that even the "RHIC" 
hydrodynamic calculations fail to reproduce azimuthally- 
integrated HBT data ^j; here, however, we simply in- 
vestigate the connection between the freeze-out geom- 
etry and oscillations of the HBT radii. Both calcula- 
tions result in a freeze-out configuration, integrated over 
PT, which is rather sharp-edged in transverse coordinate- 
space; thus, we may extract surface radii Rx and Ry to 
calculate e. 

Figure [23 shows the same quantities as plotted in Fig- 
ure [221 but extracted from these hydrodynamic calcula- 
tions. The "RHIC" source, which is geometrically ex- 
tended out-of-plane (Ry > R^, resulting in a positive 
e) generates oscillations in the transverse radii with the 
same phase as out-of-plane sources in blast-wave calcu- 
lations. For this source, the connection between e and 
the radius oscillations (Equation [221) is most robust for 
i?^ 2 aii'i least well-satisfied for -R^ 2 ' effect not ob- 
served in the blast-wave. However, our blast-wave pa- 
rameterization does not include explicit 4>p dependence 
of the temporal scale, which would affect i?o27 s.iid, to 
a lesser degree, R^^ 2- Instead of attempting a more so- 
phisticated parameterization, we simply note this fact, 
and would recommend that an experimental estimation 
of the source deformation e is probably best extracted 
from -Kg 2/-^? 0' which should be unaffected by azimuthal 
structure of the temporal scale. From the study (below) 
of parameter variations in the blast-wave, the approxi- 
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FIG. 24: (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion 1291 using several values of A'^, the maximum value of n 
taken in the summation of Equation |H| Values of the blast- 
wave parameters are for a "non-round" source, as listed in 
Table H] Arrows indicate values calculated from Equation 1321 
See text for details. 
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FIG. 25: (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion The data are not from blast- wave calculations, but 
from hydrodynamic calculations of Heinz and Kolb 1231 . for 
two different equations of state and initial conditions. Arrows 
indicate values calculated from Eauation l32l and extracted val- 
ues of freezeout edge radii; see text for details. 



mation e « 2i?2 2/i?2 ^ is good to -30%, for RHIC-type 
sources. 

Figure!^ also shows results from the "IPES" hydrody- 
namic calculation. Here, the freeze-out shape is extended 
in-plane [23, but dynamical effects are so strong in this 
extreme case, that even the sign of the transverse radius 
oscillations changes with pT- The relationships in Equa- 
tion work only at low p^, and even there only very 
approximately. According to this model, then, geometri- 
cal considerations dominate the HBT radius oscillations, 
while dynamical effects begin to dominate at much higher 
energies. 

In Figure 1261 the source geometry is azimuthally 
isotropic [Rx = Ry ^ 12.042 fm), while the flow field 
is varied from having a stronger boost in-plane (p2 > 0) 
to a stronger boost out-of-plane (p2 < 0). We notice 
again that the average HBT radius values Rf^ Q are un- 
affected by the anisotropy. The oscillations (i?^^ 2) ^-^e 
driven by flow gradients. Naively, one would expect that 



all HBT radii Ro, Ri would be smaller when the emis- 
sion angle (j)p is in the direction of the strongest boost. 
For p2 — 0.05, for example, the radii would be smaller at 
(j)p = 0, the direction of stronger transverse boost; this 
would correspond to R^ 2 < (/i = o, s, Z). For ^ — o 
and = I, this is indeed observed, at all pr- R^ 2j how- 
ever, changes sign from positive at low pT, to negative 
at high PT- This behaviour at low px is due to an ef- 
fect similar to that which led to negative proton V2 at 
low pt, even when p2 > 0. (See discussion surround- 
ing Figure IT^) In the present case, the particles with 
Px ~ are more likely to be emitted by source elements 
positioned along the y-axis, due to the strong in-plane 
boost for source elements with large spatial coordinate 
X. The homogeneity region for pT = particles is inde- 
pendent of (j}pi and has larger extent out-of-plane than 
in-plane. Thus we find the counter-intuitive result that 
Rl{(j)p = 0) > = f ) ^ RI2 > at PT « 0. We 

note that a similar argument holds for R^, except that it 
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FIG. 26; (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion l29l The anisotropy of the flow field {p2) is varied. Values 
of the other parameters are for a "round" source, as listed in 
Table H] Arrows indicate values calculated from Equation 1821 
see text for details. 



leads to the conclusion that i?^ 2 < ^-t Pt ~ 0, and so 
goes in the same direction as flow-gradient effects. It is 
only for i?^ that the two effects compete. 

Finally, comparing the scales on the right-hand panels 
of Figures 123 and it is clear that, while the 2"^-order 
coefficients are driven by both anisotropic geometry and 
flow field, a variation in geometry {Ry/R^) has a stronger 
effect on Rj^ 2 tfian a variation in p2, when these pa- 
rameters are varied by amounts which would generate a 
similar effect on elliptic flow V2 (cf Figures ITU and 
Thus, measurement of both V2 and HBT radius oscil- 
lations would allow independent determination of both 
anisotropic flow strength p2 and shape Ry/Rx- 

In Figure 071 we consider anisotropy simultaneously in 
both the source geometry and the flow field. The source 
is extended out-of-plane {Ry = 13 fm and = 11 fm), 
and the flow anisotropy (^2) varied. The 0"^-order coef- 
ficients remain unaffected by the anisotropics, while the 
i?^ 2 reflect essentially the cumulative effects of the geo- 
metric anisotropy (shown in Figure O^ and the flow field 
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FIG. 27: (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion The anisotropy of the flow field (^2) is varied. Val- 
ues of the other parameters are for a "non-round" source, 
as listed in Table Q Arrows indicate values calculated from 
Eg nation 1321 see text for details. 



anisotropy (shown in Figure I26|l , with no strong non- 
linear coupling between them. Thus, while in principle 
anisotropic flow effects may mask or dominate geomet- 
ric anisotropy 0, flow field anisotropics represent small 
perturbations on the dominant geometric effects in the 
blast-wave, using "realistic" (cf Section 4) parameters. 

In Figure|2Hlwe show the effect of increasing the trans- 
verse size of the source Ry + Rx) while keeping the 

shape (Ry/Rx) and other source parameters fixed. As 
expected, the purely-spatial transverse radius R^ (aver- 
age and oscillation amplitude) increases proportionally 
with Ry + Rx- The squared "outward" radius param- 
eter, contains both spatial components (which in- 
crease with Ry + Rl) and temporal components (which 
do not). Thus, its average value, R^q, increases almost 
proportional to Ry + Rx at low px {(3± ~ 0) and less so at 
higher px- Due to the near-cancellation of 0p-dependence 
of temporal terms, the increase in oscillation amplitudes 
^os 1 ^'^d 2 is driven mainly by the spatial terms, so 
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FIG. 28; (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion The overall transverse scale Ry + R%) is varied, 
while keeping Ry/Rx = 13/11. Values of the other parame- 
ters are for a "non-round" source, as listed in TableQ] Arrows 
indicate values calculated from Equation 1321 see text for de- 
tails. 

that -Rq, 27-^5,0 ^03,2/ ^lo display almost no sensi- 
tivity to Ry + R^. at any pT- The longitudinal radius i?^ 
is unaffected by variation in the transverse scale. 

In Figures 1291 and I3UI we vary the timescale parame- 
ters To and At, respectively. All dependence of i?^ „ on 
these parameters come directly through the dependence 
of (x^x^), which are listed explicitly in Equations 1261 
After an inspection of those equations, it is unsurprising 
that the effects of varying these parameters are similar. 
The space-time correlation coefficients which depend on 
the timescale parameters are (xi), (yi), (P), and (z^). 
According to Equation then, i?^ q and R^ 2 ^''^ unaf- 
fected by variations in tq and At. Rf = {P) is directly 
proportional to (3At^ -I- Tq), so Rf^ and i?^2 both scale 
with that quantity. Turning to the HBT radii with both 
spatial and temporal contributions, we again find that 
the <^p-dependence of the temporal terms is negligible, 
so that Rl^2/R^,Q ^-nd Rls,2/^^,o independent of the 
timescales, while R^ q displays the well-known @ sensi- 



FIG. 29: (Color online) Fourier coefficients of the (pp- 
dependence of the squared HBT radii, as calculated by Equa- 
tion 1291 The evolution duration tq is varied. Values of the 
other parameters are for a "non-round" source, as listed in 
Table H] Arrows indicate values calculated from Equation 1321 
see text for details. 



tivity to timescale. 

Thus, we find that, in the blast- wave parameterization, 
essentially all sensitivity to timescales comes through the 
(/)p-independent quantities R^o and R^q- However, it is 
important to point out that an experimental estimate 
of the freeze-out geometric anisotropy, defined in Equa- 
tion 122 from measurements of the i?^ „ would place an 
additional constraint on the evolution timescale tq. In 
particular, the large initial-state anisotropy in coordinate 
space {Ry ^ Rx 

) in a non-central collision will be reduced 
due to stronger flow in-plane than out-of-plane (p2 > 
in the present parameterization). If the source lives for 
a long time (large tq), the system may become round 
(Ry = Rx) or even in-plane extended {Ry < Rx) [23 |. 
A quantitative constraint on tq from the relationship of 
the initial to freeze-out anisotropics, however, must be 
made in the context of a realistic dynamical model and 
is beyond the scope of the blast-wave parameterization. 

In Figure 1^ the effect of variations in the (^-averaged 
("radial") flow on the HBT radius parameters is shown. 
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FIG. 30: (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion The evolution duration Ar is varied. Values of the 
other parameters are for a "non-round" source, as listed in 
Table H] Arrows indicate values calculated from Equation 1321 
see text for details. 



FIG. 31: (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion The average radial flow magnitude po is varied. Val- 
ues of the other parameters are for a "non-round" source, 
as listed in Table ^ Arrows indicate values calculated from 
Eg nation 1321 see text for details. 



As is well-known [2| , stronger flow reduces the homogene- 
ity lengths, and, indeed, almost all of the fall with 
increasing pQ. The one interesting exception is i?f 2' while 
the average scale of the longitudinal radius (i?fo) 
creases as the flow is increased, its 0p-dependence (|i?^2l) 
increases (and so, then, does \R^ 2I o\)- This is because 
there is no explicit (pp dependence in = {z^)] any 4>p 
dependence is implicit, and thus is generated by space- 
momentum correlations 44], which, in this model, arise 
solely from flow. In the no-flow limit for a boost-invariant 
source, (P), ix^), ill'^), and (z^) are all 0p-independent 
constants |44Ll47l l4a|. 

Indeed, it is for a similar reason that vanishes for 
Pt = 0, independent of model parameter in Figures 123- 
1221 At Pt = 0, symmetry demands that none of the 
spatial correlation coefficients {x^x^,) may depend on (pp. 
Hence, only HBT radii with explicit 0p-dependence may 
exhibit an oscillation in that limit. 

Finally, we note in Figure |^ that the oscillation 



strengths \R^^2\ ^^'^ 1-^8,21 ^i'^ somewhat less diminished 
(at low Pt) by increasing radial flow, than is i?^ q, which 
measures the overall spatial scale of the homogeneity re- 
gion. Here, we offer no simple insights on the interplay 
between the increasing deformation of the homogeneity 
region, and its decreasing scale, but simply note that the 
dependence of the scaled oscillation strengths on po is 
rather small, especially at low pT^ even for the unrealis- 
tic case of no average transverse flow (po = 0). 

Figure explores the effect on HBT radii of varying 
the temperature parameter T. Increasing T increases 
"thermal smearing" , reducing the flow-induced space- 
momentum correlations. It is well known that this leads 
to increased homogeneity lengths and HBT radii [s^, as 
reflected in the left panels of the Figure. In the present 
model, we find a small residual dependence of i?^ 2 ^ 
beyond the scahng of i?^ q (right panels). 

In the limit of very high temperature, all space- 
momentum correlations are eliminated, and all oscilla- 
tions of HBT radii are again due to the explicit (pp- 
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FIG. 32: (Color online) Fourier coefficients of the (pp- 
dependence of the squared HBT radii, as calculated by Equa- 
tion For T = 0.7 GeV/c, Ri o and R^q curves partially 
exceed the displayed scale. Values of the other parameters are 
for a "non-round" source, as listed in Table|l| Arrows indicate 
values calculated from Equation 1321 see text for details. 



dependence in Equations 1221 Then, we find that Equa- 
tions 1201211 hold, independent of pr- For the source of 
FigureES according to Equation^l \Rl,2/Rl,o\ = 0-083. 
FinaUy, i?^2 decreases with increasing T, and must van- 
ish at T — > cxD, when all space-momentum correlations 
are destroyed. 

Finally, we consider the effect of a finite "skin thick- 
ness" tts ■ As in the discussions surrounding Figures [7| 
and 1 171 it is appropriate to scale the flow parameters ac- 
cording to F{as) given in Eauation I2UI Moreover, it is 
clearly appropriate to scale the geometric size parameters 
Ry and Rx, as the overall scale of the source will increase 



with a. 



if these parameters remain fixed. Less clear is 

2 



the exact scaling which would keep, e.g. R^q, indepen- 
dent of Us, especially in the presence of finite flow. For 
illustrative purposes, we scale Ry and Rx also by F(as), 
so that Ry = 13.0, 12.21, and 8.84 fm, for = 0, 0.1, 
and 0.3, respectively; Ry/Rx = 13/11 in all cases. This 
scaling keeps as well as the flow gradient, the same 
in each case. 
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FIG. 33: (Color online) Emission probability contours, plot- 
ted on a linear scale, indicate emission zones for pions with 
PT ~ 0.3 GeV/c at (j}p = 0° and (jjp — 135° (indicated by 
arrows), from a blast- wave source with as — 0.3. Other pa- 
rameter values correspond to the "non-round" values listed in 
Table El 



Figure 1221 shows homogeneity regions projected onto 
the X — y plane, for a blast- wave source with as = 0.3, 
corresponding to a pseudo-Gaussian geometrical distri- 
bution in f (see Figure^. Comparing Figures ITHl and ISSl 
we observe that the lack of a "hard cut-off" in coordinate 
space in the latter case leads to less reduction in the "out" 
direction (i.e. along the direction of motion). Thus, the 
important ratio R^ o/R'i o be larger at finite pr, when 
fls — 0.3. 

We also note that the shape and size of the homogene- 
ity region itself depends much less on 0p, when = 0.3. 



The homogeneity region for pions emitted to i 



135° is 



essentially a spatially-translated (and unrotated) replica 
of that for pions emitted to (pp = 0°. Since HBT correla- 
tions are insensitive to a spatial translation of the source, 
the situation for as = 0.3 is rather similar to the situa- 
tion in which no flow is present. In the no-flow case, the 
same homogeneity region is measured at all angles, and 
all oscillations of the HBT radii arise from the explicit 
(j)p dependence in Equations 1221 43>E3; for the source in 
Figure|221 not the same region, but a (virtually) identical 
one, is being measured at all angles (pp. 

This has implications for the oscillations of the trans- 
verse radii. Focusing on the homogeneity regions for 
(j)p — 135°, we observe that while (xy) (which quanti- 
fies the "tilt" of the homogeneity region relative to the 
X and y axes) is greater when as = 0, R^s (which quan- 
tifies the tilt of the homogeneity region relative to the 
"out" and "side" directions is larger for ag = 0.3 
Therefore, we expect larger values of R^s i for larger a^. 
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FIG. 34: (Color online) Fourier coefficients of the (j)p- 
dependence of the squared HBT radii, as calculated by Equa- 
tion l29l The parameters T, tq and Ar are set to values listed 
in Table Q while the parameters po , p2, Ry and are co- 
varied with the surface diffuseness a^; see text for details. 



Figure IMI shows quantitatively the effects on the HBT 
radii, when is varied. As intended by the scahng of 
flow and size parameters with F(as), R^o remains ap- 
proximately invariant when the parameters are covaried. 
As discussed above, Rlfi/R^^o grows with as at finite 
Pt, as do the magnitudes of the oscillation strengths 
\Rl,2/Rlo\ and \Rl2/Rlo\- The average 
value of Rf remains roughly constant, while its oscilla- 
tion amplitude increases slightly with a^, due to a; — z 
and y — z correlations. 



D. Non-identical particle correlations 

Final-state interactions between pairs of non-identical 
particles (e.g. K — n) are sensitive to the space-time 
structure (size, shape, emission duration) of the emit- 
ting source, analogously to HBT correlations between 
identical particles. Most importantly, the authors of 
Ref. |63, IzS UM show that studying correlations between 
non-identical particles reveals new information about the 



average relative space-time separation between the emis- 
sion of the two particles, in the rest frame of the pair. 
This unique information may be extremely valuable to 
determine the interplay between partonic and hadronic 
effects. For example, the lower hadronic cross-section of 
kaons compared to pions may cause them to be emitted 
earlier and closer to the center of the source than pions. 
However, if most of the system evolution takes place at 
the parton level, the space-time emission pattern of pi- 
ons and kaons would be similar. This example is far from 
unique as the same kind of argument can be made for pro- 
tons recalling that pion-proton hadronic cross section are 
very large, or conversely to S and O whose hadronic cross 
section is expected to be small. Furthermore , iii terms of 
temporal effects, strangeness distillation ■Z3il23 or other 
unique physics may cause some particles to be emitted 
later than others. Preliminary analyses of the tt^ — K^, 
TT^ — , — p^ correlation functions have been re- 
ported by the STAR coUaborationjn Au-Au collisions at 
^^sWn =130 GeV and 200 GeV (Zll^- These analyses 
show that pions, kaons and protons are not emitted at 
the same average space-time point. 

The blast-wave parameterization implicitly assumes 
that the particle freeze-out conditions (temperature, flow 
profile, freeze-out time and position) are the same for all 
particle species. However, as shown earlier, the trans- 
verse momentum spectra and elliptic flow or different 
particle species do not look alike. Indeed, the relative 
contribution of the random emission (quantified by tem- 
perature) and collective expansion depends directly on 
particle masses. The same phenomena is likely to affect 
the particle freeze-out space-time emission distribution. 
In this section, we will show that collective flow effects 
implicit in the blast- wave parameterization induce a shift 
between the average freeze-out space-time point of dif- 
ferent particle species. We will then study how changing 
the blast-wave parameters affects these average freeze- 
out separations. 



1. Non-identical particle and blast-wave formalism 

Two particles interact when they are close to each 
other in phase space in the local pair rest frame. Thus, 
particle pairs may be correlated when their relative mo- 
mentum in the pair rest frame is small. For particles 
with different masses, a small relative momentum in the 
pair rest frame means that both particles have similar 
velocities in the laboratory frame and not similar mo- 
mentum. This point is particularly important to realize 
when studying correlation involving pions. Due to the 
low pion mass, the pion velocity is the same as heav- 
ier particles (e.g. kaons or protons) when its momentum 
is much lower than the particle momentum it is associ- 
ated with. For example, a proton with a momentum of 1 
GeV/c has the same velocity as a pion with a momentum 
of wO.lS GeV/c. 

As described in Ref. fl^, the spatial separation be- 
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tween particles in the pair rest frame can be projected 
along three axes, {r*^^^), along the pair transverse mo- 
mentum, (r-*j^g) , perpendicular to the pair transverse mo- 
mentum and along the beam axis. To study the 
blast- wave prediction, wc focus on the limiting case where 
the relative momentum between both particles in the pair 
rest frame is zero, which means that both particles have 
the same velocity. Thus we calculate the separation be- 
tween particle 1 and 2 in the pair rest frame {Ar*^^), 
(Ar*-^g) and {Ari^ng)^ ^ given pair transverse velocity 
Pt, azimuthal angle, (pp, and longitudinal velocity, P^: 

(Ar:„,) = - (rl,(2)) (33) 

(Ar:,,,) = - (r:,,,(2)) (34) 

(Arr„„,) = (rro„,(l))-(rr„„,(2)) (35) 

With the particle emission points defined as (x, y, z, t): 

{r*out) = lAix) cos(0p) + {y) sin(^p) - Prit)) (36) 
{r*s^de) = (y) cos(0p) - (x) sin(0p) (37) 
(ring) = 1l{{z) - pL{t)) (38) 

Recalling that {Pl) = {z/t), {rl„g{i)) = 0. Following 
the notation of Equation the variables (x) , (y) and 
(t) depend on Pt and on the particle mass (m) as: 



{x) = {x},,{m,PT) (39) 
(y) - M;^(™,/5t) (40) 
(t) = ^^' + ^" {T};;{m,PT) (41) 

To 

The average values of x and y will vary with particle 
mass and particle transverse velocity, which yields sepa- 
rations between pious and kaons, pions and protons, and 
protons and kaons in the pair rest frame shown as in 
Figure The dashed and dotted lines show the con- 
tribution of the spatial (Arout) and time (At) separa- 
tion boosted to the pair rest frame along the pair trans- 
verse momentum, respectively, i.e. the dash line shows 
{jxArout} and the dot Hue represents —{^TPrAt). The 
plain line shows (Ar*^^) = {jxiArout — PrAt)), the spa- 
tial separation in the pair rest frame. When boosting to 
the pair rest frame, time and spatial shifts in the lab- 
oratory frame add up due to their opposite signs. The 
largest (smallest) shift is obtained when the mass ratio 
between both species is the largest (smallest), i.e. for 
pion-proton (kaon-proton) pairs. 

In order to understand the behavior of Ar*„j , it is in- 
structive to investigate how the average emission points 
in the laboratory frame {{rout)) of pions and kaons be- 
have in various conditions as shown in Figure 1361 This 
figure shows the average emission points of pions and 
kaons in four different configurations: (1) The thin plain 
line shows the flow profile, (Equation^, which sets the 
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FIG. 35: Separation between the average emission point of 
pions, kaons and protons in the pair rest frame along the pair 
transverse velocity. The values of the blast-wave parameters 
are for a round source, as listed in Table |T] Dot line: time 
shift in the laboratory boosted to the pair rest frame, dash 
line: spatial shift in the laboratory boosted to the pair rest 
frame, plain line: total spatial shift in the pair rest frame. 

emission point of particles when T = GeV. (2) the 
dot line is calculated assuming an infinite system, i.e. 
n{r,(j)s) = 1 and with T = 0.1 GeV. (3) the dash line 
corresponds to a finite system as in Equation ^ with 
fls = 0, and (for illustration) an extremely low tempera- 
ture, T = 0.01 GcV. (4) the thick plain line is calculated 
using the standard parameters used in Figure|2ni Only in 
case (1) are pions and kaons emitted exactly at the same 
point. Since T = GeV, all particles are always emitted 
at the same point as set by the flow profile. In case (2), 
particle emission points spread around the average emis- 
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FIG. 36: Average emission points of pions (left) and kaons 
(right) as a function of the pion and kaon velocity and mo- 
mentum. All curve are calculated with po = 0.9 and pa =Q- 
The parameters r and /St are irrelevant. Thin plain line: 
Rx = Ry = 13 fm, T=0 GeV. Thick plain line: infinite sys- 
tem (i.e. ^{r,4>s) = 1 ), T = 0.1 GeV. Dash line: R^ = Ry 
= 13 fm, T=0.01 GeV. Dot line: R^ = Ry =13 fm, T = 0.1 
GeV. 



sion defined by the flow profile. At small 77 there are as 
many particles emitted at large r, i.e. large p as at small 
r. But when \-q\ is large, the term e,-^T/Tcoshp(r4,)cosh(n) 

in Equation ^| favors small p, i.e. small r. Thus, the 
average emission point {rout) is smaller than the one de- 
fined by the flow proflle. In case (3), the average emission 
points of pions and kaons follow closely the flow profile 
when the particle rapidity is significantly smaller than 
Po. Close or beyond poi ^ certain fraction of the par- 
ticle emission function is truncated due to the system 
finite size. The particle emission points are not allowed 
to spread beyond the system boundary, hence breaking 
the balance between particle emitted at small radii and 
particle emitted at large radii. Thus, the particle aver- 
age emission radius is smaller than the emission radius 
given by the flow profile. Because the temperature is 
rather small, the average emission points of both species 
converge rapidly toward the radius of the system. The 
larger temperature in case (4) makes the average emission 
radii converge slowly toward the system limit. This phe- 
nomenon, with the addition of the phenomenon described 
in case (2) reduces very significantly the average particle 
emission radius compared to the flow profile limit. 

The effect of temperature depends on particle mass and 
momentum. Random smearing is maximal for particles 
with low mass and momentum such as the low pT pions 
that are usually associated with kaons or protons in non- 
identical particle correlation analyses. This effect is illus- 
trated in Figure 123 It shows the probability of emitting 
pions at a given point in the transverse plane for two dif- 
ferent pion momenta ~p = (0.15, 0, 0) and 'p = (0.3, 0, 0). 
Kaon and proton momenta are calculated so that they 
have the same velocities as pions. The region of the sys- 
tem that emits particles of a given momentum shrinks 
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FIG. 37: Distribution of the emission points in the transverse 
plane, for different particle species emitted at the same veloc- 
ity, (3y = Q and = 0.907 on the left-hand side and Px = 
0.974 on the right hand side. Top panels are for pions, middle 
for kaons and bottom for protons. The same logarithmic color 
scale is used for all six panels. The blast- wave parameters are 
the same as in Figure I^Hl except that as = 0.01 The dashed 
lines show {x). 



and moves towards the edge of the system as the parti- 
cle mass, and/or momentum increases. The magnitude 
of the inward radius shift depends on the fraction of the 
source distribution that is truncated due to the system 
finite size. Thus, the inward shift of the average emis- 
sion radius scales with the source size. This effect yields 
the systematic shift between the average emission points 
of pions, kaons and protons as shown in Figure since 
pion source size is the largest and proton the smallest. 
Light particles are emitted the closer to the center of the 
source than heavier ones. 

In addition to a spatial separation, the blast-wave pa- 
rameterization induces a time shift between different par- 
ticle species emitted at the same velocity as shown by 
the dot line in Figure 1351 Due to random motion, the 
space-time rapidity (77) spreads around the momentum 
rapidity (Y). Because the relationship t = TC0sh(77) is 
positive definite, the larger dispersion of rj for pions than 
for kaons or protons leads to a delay of the emission of pi- 
ons with respect to kaons or protons. Protons are emitted 
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first, then kaons and then pions. The spatial and time 
shifts have opposite signs in the laboratory frame but 
they sum up when boosting to the pair rest frame. The 
plain line in Figure 1^ shows the added contributions of 
both shifts. The pion-kaon and pion-proton separation 
in the pair rest frame ranges from 5 to 15 fm, while the 
separation between kaons and protons is on the order of 
2-4 fm. These shifts are large enough to be measured. 

The curves on Figure have been obtained by set- 
ting the blast-wave parameters to arbitrary values. We 
will now investigate how changing these parameters af- 
fects the shift between pions and kaons. The results ob- 
tained studying pion-kaon separation can be easily ex- 
trapolated to pion-proton and kaon-proton separations. 
Since experimental analyses of non-identical two-particle 
correlations performed to date do not investigate the az- 
imuthal dependences with respect to the reaction plane, 
we will focus our study on central collisions. We will then 
show that the shift between the average emission points 
of different particle species oscillates with respect to the 
reaction plane without investigating the effect of varying 
the parameters in detail. 



2. Non-identical particle correlations in central collisions 

In central collisions, azimuthal symmetry implies that 
the particle emission pattern depends only on the relative 
angle between the position and momentum (pp—cjia- Thus, 
setting for convenience 0p = 0, yields: 

{rluS)) = {it{x,-Pt^U)) (42) 
{rUdei^)) = (43) 

Furthermore, we consider emission from an 
azimuthally-isotropic source {R^ = Ry and p2 = 0) only. 
Hence, the only quantity of interest is (Ar*„j). 

Figure I3SI shows the dependence of the average separa- 
tion between tt and K in the pair rest frame as a function 
of the pion and kaon momentum using several values of 
N, the maximum value of n taken in the summation of 
Equation |S1 Pions and kaons are of course treated as 
bosons. The difference between the average separation 
calculated either by using a Boltzman function (N=l) or 
by approximating the Bose-Einstein distribution at the 
/^th order is smaller than 0.5 fm. The maximum rela- 
tive difference is on the order of 8% at small transverse 
momentum. The Bose-Einstein distribution already con- 
verges when N=2. 

Figure 1211 shows the dependence of the spatial shift be- 
tween pions and kaons in the pair rest frame as a function 
of the pair velocity for different temperature. The shift 
increases as the temperature increases all the way from 
0.01 GeV to 0.12 GeV. When the temperature is very 
low (e.g. 0.01 GeV), pion and kaon emission patterns are 
dominated by space-momentum correlation independent 
of particle masses. In the limit of zero temperature, there 
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FIG. 38: Average separation between tt and K in the pair 
rest frame as a function of the pion and kaon momentum 
using several values of A*', the maximum value of n taken 
in the summation of Equation |H| Values of the blast- wave 
parameters are for a round source, as listed in Table Q 
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FIG. 39: Average separation between tt and K in the pair 
rest frame as a function of the pion and kaon momentum 
for different temperature. The values of the other blast-wave 
parameters are for a round source, as listed in Table |T| 



is a unique correspondence between particle velocity and 
emission point. In that case since we consider particles 
with the same velocity, pions and kaons are emitted from 
the same point. On the other hand when the temperature 
is non zero, a fraction of the pions and kaons that would 
be emitted in a infinite system are truncated, which shifts 
their average emission points inward. Because the pion 
source size is significantly larger than the kaon source size 
due to the pion lower mass and momentum, the pion av- 
erage emission point is more shifted inward than kaon's, 
which is illustrated in Figure BUI This figure shows the 
contribution of the spatial shift in the average separation 
between the pion and kaon average emission points in the 
pair rest frame. When the temperature is low (0.01 GeV), 
both spatial and time separation are small as shown by 
the dash lines in figures HUI and 1411 As the temperature 
increases, the pion emission time increases faster than 
the kaon emission time; the higher the temperature the 
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FIG. 40: Contribution of the spatial shift to the average 
separaration between pions and kaons in the pair rest frame 
as a function of the pion and kaon momentum for different 
temperature. The values of the other blast-wave parameters 
are for a round source, as listed in Table Q 



FIG. 41: Contribution of the time shift to the average sep- 
araration between pions and kaons in the pair rest frame as 
a function of the pion and kaon momentum for different tem- 
perature. The values of the other blast-wave parameters are 
for a round source, as listed in Table Q 



larger the shift between pion and kaon average emission 
time (after boosting into the pair rest frame). On the 
other hand the spatial shift varies httle within the tem- 
perature range expected in relativistics heavy-ion colli- 
sions (0.08-0.12 GeV). At temperature above 10 MeV/c, 
a fraction of the pion and kaon sources is truncated even 
at low particle velocity. The fraction of the source that is 
truncated, which leads to an inward shift of the average 
emission radius, varies with transverse momentum and 
particle mass but it is relatively insensitive to tempera- 
ture variation between 0.08 and 0.12 GeV. It is interest- 
ing to note that the average spatial separation between 
pion and kaon in the laboratory frame actually decreases 
as the particle velocity rises above 0.8c. However, this 
decrease is not visible in Figure 1^ and QUIbecayse the 
7t factor applied when boosting to the pair rest frame 
rises faster with velocity than the separation between pi- 
ons and kaons in the laboratory frame decreases. 

Changing the maximum flow rapidity (po) also affects 
the separation between pions and kaons in the pair rest 
frame as shown in Figure 1421 When po = 0, pions and 
kaons are emitted from the same space point as shown 
on FigurellSland only the time shift remains (Figure^)). 
Indeed, the time shift depends weakly on pg. Figure El 
shows that the contribution of the time shift to the sep- 
araration between pions and kaons in the pair rest frame 
reaches a plateau when the pion momentum reaches 0.25 
GeV/c. The magnitude of the time shift starts decreasing 
in the laboratory frame upon reaching this pion momen- 
tum but it is compensated by an increase of the boost fac- 
tor 7t/?t- On the other hand, when po is large enough, 
a significant spatial separation appears in the pair rest 
frame, which is sensitive to the value of po- The main ef- 
fect of increasing the flow strength is to decrease the pion 
(and kaon) transverse source size, as shown in Figure !^ 
In the laboratory frame, the spatial shift between pion 
and kaon average emission point switches from decreas- 



|i,(7t) = p^(K) = p,(pair) 



E 



^-10 



-15 



) 0.2 0.4 


0.6 


0.8 
















— Po = 

_ Po = 

Po = 

Po = 

1 




0.6 
0.9 
1.2 


< 1 . 1 < 1 < 1 . 


- 0.05 


0.1 


0.15 0.2 0.25 0.3 0.35 0.4 
p,(7i) (GeV/c) 

. < 1 < < < 1 . < < 1 < < < 1 < < . 1 



I , , , I , , I I , I , I , , , I I , , I , , , I , , I L 

0.2 0.4 0.6 0.8 1 1.2 1.4 



p,(K) (GeV/c) 

FIG. 42: Average shift between the pion and kaon emission 
point in the pair rest frame as a function of the pion and kaon 
momentum for different flow rapidity (po) The values of the 
other blast- wave parameters are for a round source, as listed 
in Table El 

ing to increasing as the particle velocity increase. The 
value of the velocity where this switch occurs, depends on 
the value of po- The increase of the spatial shift between 
pions and kaons arises from two effects: the first is the 
expected shift of the average emission point of both par- 
ticles due to the flow profile; the second is that the kaon 
source size drops more rapidly than the pion source size. 
Then, above a velocity that depends on po, the fraction 
of pion that would be emitted beyond the system limit 
starts dropping faster than the corresponding fraction of 
kaons, thus the separation between the average emission 
points of pions and kaons decreases. However, this turn 
over, which takes place in the laboratory frame is not 
directly visible in Figure US! because the boost factor jt 
compensates it. 

Figure 031 shows the sensitivity of the average separa- 
tion between pion and kaon emission point in the pair 
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FIG. 43: Contribution of the spatial shift to the average 
separaration between pions and kaons in the pair rest frame as 
a function of the pion and kaon momentum for different flow 
rapidity (po) The values of the other blast- wave parameters 
are for a round source, as listed in Table Q 



FIG. 45: Average shift between the pion and kaon emission 
point in the pair rest frame as a function of the pion and kaon 
momentum for different system radius {Rx — Ry). The values 
of the other blast- wave parameters are for a round source, as 
hsted in Table |T| 
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FIG. 44: Contribution of the time shift to the average sep- 
araration between pions and kaons in the pair rest frame as 
a function of the pion and kaon momentum for different flow 
rapidity (po) The values of the other blast- wave parameters 
are for a round source, as listed in Table Ql 



FIG. 46: Average shift between the pion and kaon emission 
point in the pair rest frame as a function of the pion and 
kaon momentum for different system proper life time (tq). 
The values of the other blast-wave parameters are for a round 
source, as listed in Table Q 



rest frame to varying the system radius. This spatial 
separation scales directly with the system radius because 
it does not modify the fraction of pions or kaons that are 
truncated due to the system finite size. 

Like the system radius, the proper life time tq acts as 
a scale driving the shift between pion and kaon emission 
time. Figure EH shows the effect of varying tq on the 
separation between pion and kaon in the pair rest frame. 
The effect of varying tq is small because the contribution 
of the time shift to the separation in the pair rest frame is 
significantly smaller than the contribution of the spatial 
separation. 

Figure El shows the effect of varying Os on the separa- 
tion between pion and kaon in the pair rest frame. Un- 
like in Figure 15^ po '^^-s kept constant. Indeed varying 
Po significantly affects the average emission radii, which 
hides the effect of changing at low velocity. However, 



when the velocity is larger than tanh(po) = 0.716c, the 
amount of boost and space-momentum correlation than 
particle acquire depends on a^; the larger Os the larger 
the boost. Thus, when the pair velocity is lower than 
0.716c, the separation decreases with increasing Os be- 
cause the fraction of truncated particles decrease. When 
the pair velocity is larger than 0.716c, increasing has 
the same consequence as increasing po- 



3. Non-identical particle correlations in non-central 
collisions 

As reported in ref [tJ, the average space-time sepa- 
ration between different particle species may depend on 
the particle emission angle with respect to the reaction 
plane. The effect in the blast-wave parameterization is 
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FIG. 47: Average shift between the pion and kaon emission 
point in the pair rest frame as a function of the pion and 
kaon momentum for surface smoothness (fls). The values of 
the other blast-wave parameters are for a round source, as 
hsted in Table Q 

shown in Figure 0S| The blast- wave parameters are the 
same as in Figure 1351 with the exception of and p2 
which are varied. Clear oscillations of {r*^f{'K)~r*^^[K)) 
are found when Rx is set to 11 fm, i.e. when the source 
is out-of-plane extended {Ry = 13 fm). Small oscillations 
of (?'*j^g(7r) — rl-^^{K)) appear as well. The oscillations 
in both directions are on the order of 1-2 fm, which may 
be measurable. On the other hand, keeping the source 
cylindrical but setting the flow modulation parameter p2 
to 0.05 yields very small oscillations, which will be very 
challenging to probe experimentally. Thus, non-identical 
two particle correlation analyses with respect to the re- 
action plane, as pion HBT may provide a good handle 
on the source shape but not on the flow modulation. 



E. Summary of effects of parameters 

Above, we have explored the sensitivity of several ex- 
perimental observables on various freezeout parameters. 
Here, we summarize in an orthogonal manner- describ- 
ing briefly the main observable effects due to an increase 
in each parameter, if the other parameters are left fixed. 

The temperature parameter T quantifies the 
randomly-oriented kinetic energy component of the 
freezeout scenario. Increasing this energy component 
leads to decreased slopes of pr spectra, especially for 
light-mass particles. Since random motion destroys 
space-momentum correlations, increasing T reduces 
measured elliptic flow (U2) and increases homogeneity 
scales, i.e. R^q. On the other hand, for the range of 
values considered, increasing temperature increases the 
average separation in the pair rest frame between pions 
and heavier particles, e.g. - r*„t(7r)). 

The (/)-averaged transverse flow strength is quanti- 
fied in this model by po. Increasing this directed en- 
ergy component decreases slopes of pT spectra, especially 
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FIG. 48: Average separation between pions and kaons in the 
pair rest frame as a function of their momentum azimuthal 
angle with respect to the reaction plane. The four different 
lines are calculated for 4 different kaon momenta and the cor- 
responding 4 different pion momenta required so that pions 
and kaons have the same velocity. The values of the blast- 
wave parameters used in the left panel are for a "non-round" 
source, as listed in Table HI. The same parameter were used 
in the middle panel except that p2 = 0. The round source 
parameters are used in the right panel except that p2 = 0.05. 
The separation along the out and side axes are shown on the 
upper and lower panel respectively. 



for heavier particles. Increasing po leads to increasing 
space-momentum correlations, which reduce i?^ g and 
{r*.ut{K) — rl^f{Tr)), and at constant p2 reduces V2 at 
high rriT. 

The parameter quantifies the "surface diffuseness" 
of the spatial density profile. Taking care to keep the 
average transverse flow fixed, variation in has little ef- 
fect on purely momentum-space observables: pt spectra 
and V2{m,pT)- On the other hand, going from a "box 
profile" {qs = 0) to a pseudo-Gaussian profile (a^ = 0.3) 
increases the "out-to-side" ratio RI^/R^.q at higher px, 
as the homogeneity region is not constrained by hard ge- 
ometric "emission boundaries." Furthermore, for a fixed 
average transverse flow (and flow gradient), increasing Os 
leads to stronger oscillations in the HBT radius parame- 
ters. Increasing Ug also increases the {rl^^iK) — r*„((7r)) 
when the temperature and particle velocity are such that 
the pion source size is significantly larger than the kaon's. 

Considering variations in the geometric transverse 



29 



scale of the source y i?^ + Ry, while keeping all else (in- 
cluding Ry/Rx) constant, we conclude that the sensitiv- 
ity is in the HBT radii and the average separation be- 
tween different particle species. In HBT, the sensitivity 
to this parameter is through R^q and R^q. The spa- 
tial shifts that contribute to the separation (r*„j(if) — 
r*„t(7r)) in the pair rest frame scales directly with this 
parameter. Momentum-space observables such as spec- 
tral shapes or V2 are unaffected. 

We turn now to variations in the timescale parameters- 
the evolution duration (tq) and emission duration (Ar). 
The only significant variations are in Rf q and R^ g , and 
these depend differently on tq and Ar, allowing indepen- 
dent probes of the two parameters. The time shifts that 
contribute to the separation {r*^^{K) — r*„f(7r)) in the 
pair rest frame scales directly with these parameters (tq 
and/or Ar). 

Non-central collisions may exhibit azimuthally 
anisotropic geometry and/or dynamics in the freeze-out 
configuration. These two effects were studied by varying 
Ry/Rx, and p2, respectively. Azimuthally- integrated 
spectra were negligibly affected, as were the 0p-averaged 
values of HBT radu (i?^_o)- Observables designed to 
probe anisotropy, V2 and Rf^ 2 of course showed strong 
sensitivity to geometric or dynamical anisotropy, in- 
creasing strongly in magnitude as Ry/ R^ deviates from 
unity and/or p2 from zero. 



IV. COMPARISON WITH EXISTING DATA 

In the previous section, we reviewed the various exper- 
imental observables that the blast-wave parameterization 
may be able to reproduce. The effect of changing blast- 
wave parameters has been examined. In this section, we 
investigate how well those parameters can be constrained 
by the experimental data. 

Before comparing the blast-wave parameterization 
with data, we highlight its limitations and possibly not 
valid assumptions: 

• Freeze-out from a thermalized system at a constant 
temperature. In the blast wave framework, freeze- 
out is assumed to take place at the same temper- 
ature for every particle. In principle, this condi- 
tion may be relaxed by fitting separately the data 
for each type of particle. (We do not attempt this 
here.) Perhaps a more important assumption, how- 
ever, is made in discussing "temperature" at all; 
it requires system thermalization. If particles are 
emitted in a non-thermal way our framework may 
fail to describe the data. 

• Longitudinal boost invariance. The blast wave pa- 
rameterization relies on longitudinal boost invari- 
ance. Elliptic flow data {v2 ) published by the PHO- 
BOS collaboration clearly show that longitudinal 
boost invariance is broken for rapidity larger than 



1 [TTJ. Thus the blast wave parameterization can 
only be applied within \Y\ < 1. This assumption 
will be validated or ruled out by high precision 
data. 

• No resonances. The observables that constrain the 
blast wave are calculated for each particle species 
using their own masses, ignoring the fact that some 
particle may be decay products of other particles. 
In fact a large fraction of pions and protons pro- 
duced at RHIC energy originate from resonance 
decays. As mentioned in section 3, pions from lo 
meson decay distort the pion source space-time dis- 
tribution but this effect is limited because the frac- 
tion of pions coming from lo is on the order of 10%. 
On the other hand, pion spectra and to a lesser ex- 
tent proton spectra may be significantly distorted 
by resonance feed-down. Within the blast wave 
framework however, it was found in |33,[z3 that ac- 
counting for resonances affects the Blast Wave pa- 
rameters extracted fit to spectra by no more than 
20%. While 20% is not negligeable, we feel that 
the gross features of the parameterization are pre- 
served. On the other hand, an investigation of the 
effect of resonance feed-down on V2 remains to be 
performed. 

We focus on data from Au-Au collisions at ^snn = 
130 GeV pubhshed by the PHENIX and STAR collabo- 
rations at RHIC. The measurements used are listed on 
top of Table IhI To avoid uncertainties between different 
experiments and to optimize the overlap between central- 
ity bins, we do not use the full set of data published at 
RHIC. Three centrality bins are available. The approxi- 
mate overlap between the centrality bins for the different 
measures may affect the quality of the fits. 

The blast-wave parameters are constrained by simul- 
taneously fitting transverse momentum spectra, trans- 
verse momentum azimuthal anisotropy (^2), and pion 
HBT radii. Bose-Einstein and Fermi-Dirac distribution 
are used when appropriate using the approximation to 
the 4*/i order as in Equation|Sl i.e. the value of n = 4 in 
the summation of Equation |S1 The quality of the fit to 
each measure is shown by the P^r data points listed in 
Table ^ The number of degrees of freedom is not shown 
because the number of parameters cannot be subtracted 
from the number of data points independently for each 
measure. In addition to the 7 blast-wave parameters, 
there are 8 scaling parameters used to normalize 7r+ , tt~ , 
, K~ , p, p. A, and A spectra. These scaling param- 
eters (A) are not minimized by the fitting function but 
calculated directly. The for a given spectra is given 
by: 

(44) 

I 

di is the measured invariant yield in the bin i. is 
the experimental error on d^, and Ci is the calculated 
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Central 


Mid-central 


Peripheral 


data 


TT, K, p spectra 


0-5% 


15-30% 


60-92% 


A spectra [80j 


0-5% 


20-35% 


35-75% 


pion radii 


0-12% 


12-32% 


32-72% 


Elliptic flow [38] 


0-11% 


11-45% 


45-85% 


X^/(# data points) 


TT^ & 7r~ spectra 


7.2/10 


26.5/10 


13.0/9 


/sT"*" & K~ spectra 


24.2/22 


21.4/22 


10.1/10 


p &ip spectra 


10.6/18 


23.2/18 


28.0/12 


A & A spectra 


9.5/16 


12.8/16 


11.0/16 


TT 172 


14.6/12 


29.8/12 


5.2/12 


P V2 


1.6/3 


9.2/6 


0.8/3 


Tout 


1.9/6 


0.4/2 


0.4/2 


^ T' side 


2.7/6 


0.07/2 


0.06/2 


7^ Tlong 


5.3/6 


0.003/2 


0.1/2 


Total 


77.6/99 


107.7/90 


68.7/68 


parameters 


T (MeV) 


106 ± 3 


107 ± 2 


100 ± 5 


Po 


0.89 ± 0.02 


0.85 ± 0.01 


0.79 ± 0.02 




0.52 ± 0.01 


0.50 ± 0.01 


0.47 ± 0.01 


P2 


0.060 ± 0.008 


0.058 ± 0.005 


0.05 ± 0.01 




13.2 ± 0.3 


10.4 ± 0.4 


8.00 ± 0.4 




13.0 ± 0.3 


11.8 ± 0.4 


10.1 ± 0.4 


T(/m/c) 


9.2 ± 0.4 


7.7 ± 0.9 


6.5 ± 0.6 


At(/m/c) 


0.003 ± 1.3 


0.06 ± 1.3 


0.6 ± 1.8 



TABLE II: Upper section: data used in the fit. Middle sec- 
tion: number of X^/ data points for each measure. Lower 
section: best fit parameters. Note that (/3t) is not a fit pa- 
rameter, but it is calculated from po- 
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FIG. 49: Comparison of the data with the blast-wave calcu- 
lations performed with the best fit parameters in three cen- 
trality bins. The closed circles are central data, the open cir- 
cles are mid-central data and the crosses are peripheral data. 
The plain lines show the blast wave calculation within the fit 
range while the dash lines show the extrapolation over the 
whole range. 



invariant yield. The minimum occurs when dx'^/dA 
= 0, hence: 

To avoid the region where resonance feed-down may 
be large and where the Boltzmann approximation may 
not be valid, fits to transverse momentum spectra are 
restricted to pions with pr > 0.4 GcV/c. Wc further 
restrict the fit range to pt < 1 GeV/c for pions, px < 1-5 
GeV/c for kaons, pr < 2 GeV/c for protons and lambdas, 
because hard processes may significantly contribute to 
the particle yield above these momenta zS\. 

Given the large number of data points from spectra, 
the temperature (T) and the maximum flow rapidity (po) 
parameters are constrained by transverse mass spectra. 
A comparison of data and the blast- wave calculation with 
the best fit parameters is shown in Figure 021 In semi- 
logarithmic plots, the data and the model fits seem to 
agree very well. The x^/(# data points) in Table UTI 
show, however, that the fit quality varies between particle 
species and centrality bins. The correlation between the 



flow rapidity and temperature is shown in Figure lSUI The 

extracted temperatures and flow velocities are consistent 
with similar studies reported in Rcf. 29, 30, 32, 81]. The 
differences between these analyses are due to the use 
of different data sets, different fit ranges, and different 
flow and spatial profiles. The study reported in Ref. [s^l 
shows that accounting for resonance feed-down leads to 
a temperature 30% higher and a transverse velocity 15% 
lower than in our study. 

The fit of V2 mainly constrains the flow rapidity modu- 
lation parameter, p2, and the ratio between the radius in- 
plane and out-of-plane, Rx/Ry A comparison between 
the data and the blast-wave calculation with the best 
fit parameters is shown in Figure 1511 The correlation 
between p2 and Rx/Ry is shown in Figure 122 The fit 
performed in this paper differs from the fit reported to- 
gether with the data in Ref. 38J. There are three main 
differences: (i) in Ref. 0], the temperature and flow ve- 
locity are free parameters while in this paper they are 
mostly determined by spectra shapes, (ii) a filled vol- 
ume and a flow profile are used instead of a shell, (iii) 
the fits are performed on three centrality bins instead of 
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FIG. 50: Number of a contour of the maximum flow rapidity 
(po) with respect to the temperature (T). The shadows show 
the 1, 2 and 3 a contours from darkest to hghtest. 
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FIG. 51: Comparison of the «2 data with the blast-wave cal- 
culations performed with the best fit parameters in three cen- 
trality bins. The closed circles are central data, the open cir- 
cles are mid-central data and the crosses are peripheral data. 
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FIG. 52: Number of a contours of the ratio of the in-plane 
radius over the out-of-plane radius (Rx/Ry) with respect to 
the flow rapidity oscillation parameter (p2). The colors show 
the 1, 2 and 3 a contours from darkest to lightest. 
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the minimum bias sample. Thus, a direct comparison 
between the fit parameters extracted in this paper and 
those reported in Ref. [s^ is difficult. The freeze-out 
spatial anisotropy probed by the ratio Rx/Ry, increases 
when increasing the initial state spatial anisotropy, i.e. 
when going from central to peripheral events. The flow 
modulation decreases from central to peripheral events, 
following the same trend as po- 

The pion source radii allow us to estimate the pa- 
rameters Rx (or T and At. As discussed earlier, 
the ratio Rx/Ry is constrained by the fit to V2. In the 
mid-peripheral and peripheral bins, the pion source size 



FIG. 53: Comparison of the pion source data measured by 
the STAR (circle) and the PHENIX (box) collaborations with 
the blast-wave calculations. Only the STAR data were used 
to constrain the blast-wave parameters. The closed circles are 
the central data, the open circles are the mid-central data and 
the crosses are peripheral data. 



has been measured at only one transverse momentum, 
which weakly constrain the parameters. The compar- 
ison between the data and the blast-wave calculation 
performed with the best fit parameters is shown in Fig- 
ure l53l This figure includes the pion source radii reported 
by the PHENIX collaboration |51| although they weren't 
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used in the fit. These radii were measured over the 30% 
most central events, which doesn't overlap well with the 
STAR collaboration centrality bins (0-12%, 12-32% and 
32-72%). Qualitatively the PHENIX data points agree 
with the blast-wave fit as they oscillate between the cen- 
tral and mid-peripheral lines. The obtained from the 
fit is small, because the systematic errors were added 
in quadrature to the statistical errors. The system ra- 
dius increases from peripheral events to central events 
following the growth of the system initial state. The final 
system radius is roughly twice the initial system radius, 
which is a clear evidence of the system expansion. The 
proper time (r) is surprisingly small. Typically, hydrody- 
namical calculations reach kinetic freeze-out in 15 fm/c, 
not 9 fm/c '82']. The hydrodynamical calculations may 
be over-predicting the system life-time or the assump- 
tion underlying the extraction of r may not be valid. 
The assumption of a boost invariant expansion along the 
beam axis from a narrow (less than 1 fm thick) pancake 
may be too simplistic. Indeed, it has been shown that ac- 
counting for the fluid viscosity breaks the boost invariant 
flow |25l |. Furthermore, three dimensional hydrodynami- 
cal calculations reproduce the measured pseudo-rapidity 
distribution only when the initial width of the system 
along the beam axis extends over several units of spa- 
tial pseudo-rapidity "ss]. Such initial conditions lead to 
a longitudinal expansion that also breaks boost invariant 
scaling. In these calculations, even though the system 
persists for 15 fm/c, the Riong radii is well reproduced. 
Thus, the absolute value of proper life-time, r, extracted 
within the blast-wave parameterization may be question- 
able. Along the same lines, the very short emission du- 
ration is surprising. Indeed, a scenario where particles 
are emitted in a flash would appear rather unlikely. (On 
the other hand, processes such as supercooling |8J] have 
been discussed, and may lead to such a scenario.) Re- 
cent experimental data show that the value of the 
Rout radius may have been underestimated by the exper- 
iments. Thus, since the emission duration parameters At 
strongly depends on Rout, it is not appropriate to draw 
definite conclusions until the experimental issues are set- 
tled. 

We have shown that tt, K, p, A transverse momentum 
spectra, tt, p elliptic flow, and pion source radii measured 
in Au-Au collisions at y^sjvw — 130 GeV are well repro- 
duced by the blast-wave parameterization. With the ex- 
ception of the time-scale parameters (system proper life 
time and emission duration), which fall short of any re- 
alistic model calculations, the fit parameters are within 
expectations. This issue may be resolved when new data 
on pion source radii become available. Based on the pub- 
lished data in Au-Au collisions at ^^snn = 130 GeV, we 
conclude that the blast-wave parameterization provides 
a good description of the system freeze-out stage. 

This conclusion will be tested in the future, using the 
sensitivity of the blast-wave parameterization to observ- 
ables that have been presented in the previous section: 
the oscillation of the pion radii with respect to the re- 
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FIG. 54; Oscillations of the pion source radii obtained with 
the best fit parameters in the three centrality bins. 



action plane, and the space-time separation between the 
emission points of different particle species. Using the 
parameters extracted from the fits we have calculated 
the corresponding oscillation of the pion radii with re- 
spect to the reaction plane (Figure 154(1 and the sepa- 
rations between the average space-time emission point 
of pious kaons and protons (Figure I55|l . We have also 
calculated the kaon source radii as shown in Figure 1561 
since they inay become available from the RHIC exper- 
iments |85l I8(tI | . The blast- wave parameterization faces 
the challenge of simultaneously reproducing a large vari- 
ety of observables that will be measured in Au-Au col- 
lisions at 200 GeV, namely: (1) pion, kaon, proton and 
A transverse momentum spectra, (2) the elliptic flow of 
many particle species, (3) the pion source radii including 
the oscillations with respect to the reaction plane, (4) the 
kaon source radii, (5) the space-time separation between 
pion, kaon and proton sources. If a satisfactory agree- 
ment between these various measurements is achieved, 
it would provide evidence of a collective expansion that 
would be very challenging to avoid. Furthermore, when 
high quality data becomes available that allows for sepa- 
rate fits to each particle species, the blast-wave parame- 
terization may prove to be an important and handy tool 
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FIG. 55: Space-time separation between pions and kaons 
(left), pions and protons (middle) and kaons and protons 
(right) calculated with the best fit parameters in the three 
cciitrality bins. Plain line: parameters from fit to central 
data. Dash line: parameters from mid-pcriphcral data. Dot- 
ted line: parameters from fit to peripheral data. 
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FIG. 56: Kaon source radii (plain line) compare to the pion 
radii (dash line). The blast- wave calculation are performed 
with the best fit parameters in each centrality bins. 



to assess whether or not the freeze-out time depends on 
the hadronic cross-sections. It will be especially interest- 
ing to study particles with presumably low hadronic cross 
sections such as the 0, !B, f2. Since these particles may 
not be sensitive to hadronic collectivity, they may have 
a higher freeze-out temperature, they may pick up less 
collective flow, and freeze-out earlier and from a smaller 
system than tt, K, p, A. It may be possible to distinguish 
between the collective expansion at the partonic (if any) 



versus hadronic stage. 



V. SUMMARY AND OUTLOOK 

We have discussed a scenario of particle kinetic freeze- 
out in heavy ion collisions, inspired by the results of full 
hydrodynamic calculations of the collision evolution. The 
scenario, while simplified, nevertheless includes several 
important features which drive experimental observables 
of the bulk properties of the system. This includes ran- 
dom ("thermal") motion superimposed upon collective 
flow fields, an anisotropic transverse geometry, and a 
freeze-out distribution in proper time. The main assump- 
tions of the model are: longitudinal boost-invariancc; the 
same parametric source for all particle types (e.g. n,p); 
and the invariance of source parameters over the kinetic 
freeze-out process (or, cquivalently, that the model source 
parameters can represent the appropriate average of the 
time-evolving parameter). We discussed the interplay of 
the various features, such as space-momentum correla- 
tions depending on the competition between thermal and 
collective motion and the source geometry. 

The various features (e.g. random thermal motion) 
were quantified with several parameters (e.g. T) in a 
mathematical model outlined in Section 2. The general 
class of integrals which relate the model emission function 
to the various experimental observables was identified. 

The sensitivity of several observables to the underly- 
ing physics represented by the parameters was studied 
by systematically varying each parameter in Section 3. 
In addition to momentum-space observables such as the 
Pt spectra and elliptic flow, which have been discussed 
previously in the context of similar models, we focussed 
on the sensitivity of final-state correlations between non- 
identical particles and azimuthally-sensitive HBT, new 
coordinatc-spacc-scnsitivc; analysis tools now becoming 
available with the quality high statistics datasets from 
RHIC. 

The non-trivial interplay between competing physics 
effects was different for the different observables. It 
is clear that the model parameter-space can only be 
constrained by examining several observables simulta- 
neously; e.g. Pt spectra alone are insensitive to the 
source shape {Ry/Rx) and only crudely distinguish be- 
tween high transverse flow and high temperature. Fur- 
thermore, from Section 3, it is clear that experimental 
data can ofer-constrain the model parameter space. If 
consistency with px spectra and HBT data demand a 
high temperature and small transverse boost velocity, 
this has inescapable conseqiienccs for tt — K correlations 
and elliptic flow. This allows for the "break-down" of 
the parametric model, and indicates the need to consider 
new driving physics effects. 

In Section 4, the model was used as a functional form 
in a simultaneous flt to particle spectra (^^(Pt, ti)), 
elliptic flow {v2{pT, m)), and pion HBT {R{pt)) for three 
centrality classes from published RHIC results for Au-Au 
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collisions at ^/sWn = 130 GeV. While higher-pT HBT 
radii from PHENIX are only qualitatively described, a 
good fit is obtained to the pr spectra, elliptic flow, and 
\ow-pT STAR HBT radii. 

Most of the physical parameter values obtained from 
the fits, as well as their evolution with collision cen- 
trality, fall within reasonable expectations. The excep- 
tions are the evolution duration tq and emission dura- 
tion At, which are shorter than most physical models 
predict. The At discrepancy may be partially resolved 
by an improved experimental treatment of the Coulomb 
interaction in HBT analyses; this will be evaluated upon 
publication of the new results from RHIC. The tq value 
is strongly connected to the model assumption of boost- 
invariance; if this assumption is invalid at mid-rapidity 
at RHIC (e.g. if there is longitudinal acceleration in the 
early system evolution), we expect our fit values to be 
only lower limits on the evolution duration. 

Finally, Section 4 discusses expectations for kaon HBT, 
azimuthally-sensitive pion HBT, and emission space-time 
separations between non-identical particles (tt, K, p), 
based on the fit parameters to the published data. Ex- 
pectations for the azimuthally-sensitive pion HBT are 
driven by the elliptic flow data, which indicates Ry > Rx 
(see Table i.e. the freeze-out system shape is out- 
of-plane extended, qualitatively similar to the geometri- 
cal overlap configuration in the entrance channel of the 
collision. From Section [HI CI it is clear that this deter- 
mines the phases (and amplitudes) of the expected oscil- 
lations in the HBT radii- i.e. the signs and magnitudes 
of i?^ 2 ■ If azimuthally-sensitive HBT measurements con- 
firm this out-of-plane freeze-out configuration, it might 
provide further evidence of short evolution times tq, since 
it would imply that the source, which is expanding faster 
in the reaction plane than out of it, did not have suf- 
ficient time to overcome its initial geometric anisotropy 
before freeze-out. Quantifying the timescale in this way 
is outside the scope of the blast-wave model, which does 
not attempt to describe the system evolution. 

Exploration of simple parameterizations such as the 
blast- wave is driven by the desire to connect observations 
to driving underlying physics. Furthermore, a quantita- 



tive interpretation of many measurements (perhaps espe- 
cially two-particle interferometry) requires some model 
assumptions. Naturally, one prefers assumptions which 
are consistent with other measurements (e.g. elliptic 
flow) in the same system. Especially since first-principle 
model calculations have difficulty consistently explaining 
the range of observations in the soft sector of RHIC, it 
is hoped that insight may be gained from rather simple 
theory- mspirerf parameterizations such as the blast-wave. 

Inasmuch as such parameterizations reproduce obser- 
vations, there is hope that the driving physics effects have 
been approximately quantified, and one may compare 
the "underlying" parameters (e.g. timescales) with first- 
principles calculations in an effort to isolate the cause of 
the discrepancy between such calculations and observa- 
tion. On the other hand, when simplistic parameteriza- 
tions cannot, with any set of parameters, reproduce the 
main features of the data, it indicates that something 
additional is driving observations. In parameterizations, 
features may be implemented and "turned on and off," 
one at a time, testing for consistency with a range of 
experimental observations. The insight gained from such 
studies can then be fed back into more "theoretical" mod- 
els. It is hoped that new experimental data will further 
constrain- or break- the simple parameterization, gener- 
ating insights into the dynamics of heavy ion collisions 
at RHIC. 
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